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Abstract 



In the present work we shall describe and apply the techniques of the Renormalization Group - 
based in data rescaling and operator renormalizing - and of Homogenization - that substitutes, 
in a certain limit, a periodically inhomogeneous médium by a homogeneous one - for the study 
of asymptotic behavior of solutions of PDE's. 

One class of problems studied in this dissertation consists of nonlinear, diffusive differential 
equations with periodic coefficients, that supposedly model the diffusion in a heterogeneous 
médium. Our study of these problems involves arguments originated from the two techniques 
described above. This makes, in certain sense, a connection between both of them. 

Another class of problems studied here consists of nonlinear, diffusive equations with time de- 
pendent coefficients, originated from montecarlo simulation of the asymptotic behavior of the 
averaged solutions of stochastic differential equations modeling the phenomenon of two-phase 
flow in porous media. 

Using the Renormalization Group we shall classify the problems in accordance with their qual- 
itative behavior. This classification will be verified numerically. We shall also study the curve 
of phase transition that separates asymptotic regimes predominantly linear and nonlinear. The 
main results presented in this work are given by a numerical version of the Renormalization 
Group, which will be described in detail in the body of this dissertation. 
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Resumo 



Neste trabalho vamos descrever e aplicar as técnicas do Grupo de Renormalização - baseada em 
mudança de escalas e renormalização de operadores - e de Homogeneização - que substitui, num 
certo limite, um meio periodicamente não-homogêneo por um outro homogéneo - para o estudo 
do comportamento assintótico de soluções de EDP's. 

Uma classe de problemas estudados nesta dissertação consiste em equações diferenciais difusivas 
não-lineares com coeficiente de difusão periódico, que supostamente modelam a difusão em um 
meio heterogéneo. O nosso estudo desses problemas envolve argumentos originados das duas 
técnicas descritas acima, fornecendo, em certo sentido, uma conexão entre ambas. 

Outra classe de problemas estudados nesta dissertação consiste em equações difusivas não-lineares 
com coeficiente de difusão dependente do tempo, que têm sua motivação física dada pelo compor- 
tamento assintótico médio, com respeito a simulações estocásticas, do fenómeno de escoamento 
bifásico em meios porosos. 

Utilizando o grupo de renormalização vamos classificar os problemas de acordo com seu com- 
portamento qualitativo. Esse comportamento qualitativo será verificado numericamente. Vamos 
também estudar a curva de transição de fase entre os regimes assintóticos predominantemente 
linear e o não-linear. Os principais resultados apresentados neste trabalho são obtidos via uma 
versão numérica do grupo de renormalização, que será descrita em detalhes no corpo desta dis- 
sertação. 
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Capítulo 1 



Introdução 



1.1 Auto-similaridade e universalidade 



A evolução no tempo de sistemas físicos fora do estado de equilíbrio é muitas vezes bem descrita 
por soluções de equações diferenciais parciais. Frequentemente é observado que, sob circun- 
stâncias apropriadas (mas razoavelmente gerais), a evolução temporal de tais sistemas torna-se 
assintoticamente auto-similar 1 . No contexto de EDP's, tal afirmação traduz-se no fato de a 
equação em questão ter soluções que se comportam como 



Ocorre muitas vezes que tal comportamento auto-similar possui um forte viés de universalidade: 
os expoentes a e (3 e a maioria dos aspectos da função </>(•) independem do problema de valor inicial 
associado, isto é, independem dos dados iniciais (dentro de uma certa classe), dos parâmetros da 
equação ou até mesmo da sua forma. 

A análise de comportamentos assintóticos auto-similares envolve um problema em escalas múl- 
tiplas. A palavra escala se refere à escolha de uma unidade de medida, tanto para a variável 
espacial x quanto para a variável temporal t. Estas escalas são interdependentes e em geral elas 

surgem de uma determinada simetria da EDP em questão - ou parte dela. Essa simetria consiste 

1 significado de "auto-similar" está relacionado com a invariância com respeito à ação de um certo grupo de 
simetria, como será explicado logo adiante. 




quando t — > oo. 
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no fato de que existem números reais a e (3 e uma função (/)(■) tais que a equação admita uma 
solução da forma 

t*(M) = ^(£)- (1.2) 

Essas soluções são ditas auto- similares ou invariantes por mudança de escala. Uma função 
homogénea de grau —a é o exemplo mais simples de uma função auto-similar. De fato, considere 
a mudança de escala t — > L t, onde L é um fator positivo, e considere o conjunto C das funções 
/(•) : M.* + — > IR. Dado a G IR, defina, sobre C, a operação de mudança de escalas f(t) = L a f(Lt). 
Se y(t) = At~ a , t > 0, então y(t) G C e = y(t), isto é, ?/(í) é invariante pela mudança de 
escala em C. 

De forma análoga, assumindo a existência dos expoentes a e (3 em (I1.2J1 e dados L > e 
u : IR x IR^ — > IR, definimos a mudança de escalas u > v por 

v(x,t) = L a u(L p x,Lt). (1.3) 

Claramente, soluções da forma (J1.2|) são invariantes por esta mudança. Os expoentes a e (3 
são ditos expoentes críticos e a função é chamada função perfil. Em alguns casos, a equação 
completa pode não ter tal simetria mas, sob certas condições, ela ainda é "assintoticamente" 
simétrica e a solução do problema de valor inicial se comporta, para tempos longos, como 

u(x,í) w^(^). (1-4) 

O prefator A é supostamente não-nulo e pode depender da equação diferencial em questão (dos 
parâmetros e dos termos não-lineares interpretados como perturbações) e da condição inicial. 
Contudo, os expoentes críticos e a função perfil são, muitas vezes, universais, isto é, eles não de- 
pendem nem do dado inicial e nem das perturbações adicionadas à equação, desde que escolhidas 
dentro de uma classe adequada. 

Definimos classes de universalidade como conjuntos de problemas que apresentam um mesmo 
comportamento assintótico. Comportamentos da forma (jl.4|) são completamente determinados 
por: expoente a, expoente /3, função de perfil e prefator A. A classe de universalidade será 
definida em termos de a, (3 e (b. 
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Por exemplo, dizemos que duas condições iniciais estão na mesma classe de universalidade se 
as soluções dos respectivos problemas de valor inicial tiverem um comportamento assintóticoda 
forma (jl.4|) com os mesmos expoentes a e j3 e mesma função de perfil 0. 

Estendendo a noção de expoente crítico, diremos que a e (3 são expoentes críticos se u satisfizer 

u(x,t) * ^-<f> ( , (1.5) 



gi(ty \g 2 (t) 

onde g± e g<i são funções positivas tais que ln gi (t) — > a e ln g 2 (t) —> j3 quando t — > oo. A 
relação (|1.5|) é mais precisamente formulada através da existência do limite 

9i(t) u(g 2 (t) ■ , t) — ► A </>(•). (1.6) 

Essa definição engloba tanto o comportamento descrito por (jl.4|) como também outros casos de 
interesse, como, por exemplo, 



(t\nt) a xtf 3 

1.2 Problemas estudados 



Neste trabalho estaremos interessados em estudar o comportamento assintótico de soluções de 
problemas de valor inicial difusivos. Listamos abaixo os problemas de interesse. 

Problema 1 (Difusão Não-Linear com Coeficiente Periódico) 

u t = (1+ (J,g(x)) u xx + XF(u, u x , u xx ) 
u(x,0) = f(x), 

onde g é continua e periódica 2 , /j, G R e tal que 1 + f^g{x) > Vx G R e F é uma perturbação. 
A condição inicial f(x) é suave e de suporte compacto. 

Se fi = temos a chamada equação de difusão não-linear. Em particular, se F(u, u x , u xx ) = 
—uu x + perturbações de ordem superior, então a equação é conhecida como equação de Burgers 



2 Sem perda de generalidade, vamos assumir que g tem período 2n. 
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não-linear. Ainda, se a perturbação de ordem superior depender somente de u, isto é, se F for 
da forma F(u, u x , u xx ) = —uu x + g(u), então a equação é dita de reação-difusão-convecção. 

No caso geral (// 7^ 0), temos um modelo aproximado para o fenómeno da difusão em um 
meio composto e condutor, caracterizado por uma estrutura microscópica periódica, além das 
não-linearidades mencionadas acima. Esse meio é não-homogêneo, pois o coeficiente de difusão 
K(x) = (1 + (ig(x)) não é constante com respeito à variável x. 

Problema 2 (Difusão Não-Linear com Coeficiente Dependente do Tempo) 

f u t = K(t) u xx + XF(u, u x , u xx ) 
\u(x,0) = f(x) 

onde 

K(t)=t p + o(t p ), (1.8) 

p > 0, o(s) sendo uma ordem pequena de s, ou seja, lim — — = 1. A condição inicial f(x) é 
suave e de suporte compacto. 

Equações desse tipo, em que a difusividade cresce com o tempo, aparecem na descrição de certos 
escoamentos multifásicos em meios porosos. Entre estes se incluem escoamentos associados a 
processos de recuperação de hidrocarbonetos em reservatórios de petróleo, e de remediação e 
controle de poluição em aquíferos |2*5] . Nestes escoamentos, o campo de velocidades é propor- 
cional ao gradiente de pressão, segundo a lei de Darcy, e o fator de proporcionalidade é uma 
propriedade geológica do meio poroso, denominada "permeabilidade" (recíproca da resistividade 
ao escoamento). Devido à sua grande variabilidade como função da posição no meio poroso, e à 
impossibilidade de uma caracterização precisa desta variabilidade, a permeabilidade é comumente 
modelada como uma função aleatória da posição (campo estocástico) [Hj. Em consequência, as 
equações que descrevem o escoamento possuem coeficientes estocásticos e, portanto, suas soluções 
sao estocásticas. Ou seja, previsões baseadas nestas equações são inerentemente estocásticas. 

Nesse contexto, a equação acima descreve o escoamento esperado (médio, no sentido estatístico) 
na presença de heterogeneidades geológicas estocásticas [THE)]. A difusividade crescente como 
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função do tempo resulta da ocorrência dessas heterogeneidades em todas as escalas de compri- 
mento. 

Em geral, o comportamento assintótico das soluções dos problemas descritos acima é auto-similar 
- isto é, da forma (jl.4|) - ou quase auto-similar, com uma correção logarítmica - isto é, da 
forma (jl.7|) - e possui características bastante interessantes. Enfatizamos algumas abaixo. 

Os expoentes a e (3 são universais: independem da condição inicial /, da perturbação periódica 
Hg, do coeficiente K{t) e da perturbação não-linear F, desde que escolhidos dentre uma classe 
adequada (nesta classe, a e (3 dependem apenas do expoente p em (ll.Hjl ). A função de perfil 
também é universal dentro desta classe: no Problema ^ depende apenas da média harmónica do 
coeficiente de difusão, enquanto no Problema 121 depende apenas do expoente p. 

No caso linear e homogéneo (/i = = A no Problema Q A = no Problema |2J) , o prefator A 
é universal, dependendo apenas do que chamamos de massa do dado inicial /. Mais especifi- 
camente, temos A = J R f(y)dy, ou seja, A não depende de nenhum dos parâmetros da equação 
nem de detalhes em pequena escala do dado inicial /. No caso não-linear, este prefator não é 
universal: o que os teoremas ou os estudos numéricos fazem é mostrar a existência de algum A 
estritamente positivo para o qual vale o limite (jl.fij) ou, no máximo, determiná-lo. 

1.3 A técnica do grupo de renormalização 

A técnica que nos permite estudar os problemas aqui apresentados e vários outros semelhantes 
chama-se a técnica do grupo de renormalização, que daqui em diante chamaremos de RG 3 . Do 
ponto de vista de Equações Diferenciais, a técnica foi desenvolvida no final dos anos 80 e início 
dos anos 90 por Goldenfeld e colaboradores (veja fH] e referências lá citadas), em sua versão 
formal, e por Bricmont et al. [HE], em sua versão rigorosa. A versão numérica foi introduzida 
por Chen e Goldenfeld jH]. 

A aplicação formal desta técnica constitui-se, grosso modo, dos seguintes passos. Primeiro, 
3 A sigla RG vem da expressão Renormalízation Group, do inglês. 
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procuramos por uma solução da EDP que seja da forma (jl.2|) . Se a EDP contém termos não- 
lineares, nós os ignoramos ou ignoramos parte deles e consideramos apenas a parte da equação 
que nos interessa. A suposição de que a equação possui solução da forma (|1.2|) é muito forte e leva 
à inevitável redução do grau de liberdade do problema. Isso significa que reduzimos a EDP a uma 
EDO, cuja incógnita é <p : IR — > R. Se formos capazes de encontrar tal solução teremos a solução 
auto-similar desejada. O próximo passo é descrever a classe de universalidade desta solução, ou 
seja, o conjunto de dados iniciais e perturbações para os quais o comportamento assintótico é 
ditado por esta solução auto-similar. Consideramos a mudança de escalas (jl.3|) e renormalizamos 
a equação diferencial em questão, bastando para isto obter a equação satisfeita por v. A parte 
do operador diferencial que consideramos ao encontrar a solução auto-similar deve ser invariante 
por esta renormalização, ao passo que as partes restantes, às quais chamamos de perturbações, 
em geral sofrem alterações. Classificam-se então, formalmente, as perturbações como irrelevantes 
no caso em que a sua contribuição no operador diferencial tende a se anular com os sucessivos 
processos de renormalização; como marginais as perturbações que permanecem invariantes ou 
assintoticamente invariantes pela renormalização; como relevantes as que tendem a crescer neste 
procedimento. Repetimos esta mudança de escalas acompanhada da renormalização da equação 
por várias vezes e esperamos verificar que, sob certas condições, essas iterações resultem em 
alguma convergência. Por isso, a técnica do Grupo de Renormalização é considerada uma técnica 
em escalas múltiplas. 

A aplicação rigorosa do método segue uma idéia simples e elegante: ao invés de tomarmos o 
limite t — > oo de uma só vez, estabelecemos um procedimento iterativo em que o limite t — > oo 
é visto como a composição de um operador consigo mesmo um número infinito de vezes. Esse 
operador chama-se transformação do grupo de renormalização e a sua definição depende da 
equação em questão. Os pontos fixos da parte linear do operador são as soluções invariantes 
por escala. Determina-se então a bacia de atração de cada ponto fixo. Para condições iniciais 
dentro da bacia de atração, o fluxo do RG fornecerá então o comportamento para tempos longos. 
Resumindo, o procedimento se reduz aos seguintes passos: dada uma escala L > 1, integram-se 
as equações, evoluindo o dado inicial do tempo 1 ao tempo L 2 ; reescalonam- se a variável espacial 
x e a solução u por fatores de L (esses fatores dependem dos expoentes críticos); subtrai-se da 
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solução reescalonada a contribuição vinda do ponto fixo, o resto supostamente é "irrelevante" e 
espera-se que seja contraído à medida que iterarmos esse procedimento; a solução reescalonada 
é então utilizada como condição inicial para uma nova equação, similar à original, dita equação 
renormalizada; o procedimento é iterado até que o regime assintótico seja atingido. 

Um exemplo da aplicação do RG analítico é o estudo para a versão homogénea do Problema Cl 



isto é, com /j, = 0. Bricmont et al. [Zj mostraram que, para F da forma F(u,u x ,u xx ) 



u a ului„ 



e para expoentes a, b e c tais que a + 26 + 3c — 3 > 0, se a condição inicial / for suficientemente 
pequena então vale 

u(x,t) ps ~^4>* ' 1 uan do t^$> 1, 



onde A é uma constante positiva e 



0*0) = -4= exp (-%- 1 . (1.9) 



Vã^ V 4 

O sentido do símbolo ps acima é o de um limite da forma mas o resultado acima é provado na 
topologia de um certo espaço de Banach cuja formulação exata não nos interessa neste momento. 
O que chama a atenção é o fato de que, dentro de um conjunto consideravelmente amplo, o 
teorema fornece um critério extremamente simples que classifica aquelas perturbações F que não 
modificam o comportamento assintótico da solução, em relação àquele da equação do calor (caso 
em que, além de /i = 0, temos À = 0). Fazendo 

d F = a + 26 + 3c — 3, 

Bricmont et al. [Zj classificam as perturbações como irrelevantes se dp > 0, marginais se dp = 
e relevantes se dp < 0. O resultado ora mencionado será apresentado de forma precisa mais 
adiante (veja Teorema 12.71 no Capítulo |2J). Observe também que o teorema nos diz que os 
expoentes críticos a e j3 valem 1/2 e que a função perfil é a distribuição gaussiana, seja qual for 
a perturbação irrelevante adicionada à equação linear. 

No Problema El a classificação formal das perturbações é semelhante à do Problema d mas no 
lugar de dp temos 

ri F = a + 26 + 3c- (1.10) 
p+ 1 

Diremos que o termo não-linear é irrelevante se tjf > 0; marginal se t]p = 0; relevante se rjp < 0. 
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A primeira novidade deste problema é que qualquer p > faz com que as perturbações da forma 
u 3 , uu x e u xx deixem de ser marginais para se tornarem irrelevantes (repare que dp = para 
essas três perturbações, ou seja, no caso do Problema [T] essas perturbações são marginais). 

Usando as mesmas ideias introduzidas originalmente por Bricmont et al. [7j, Braga et al. [3] 
provaram que as soluções do Problema El linear, isto é, com A = 0, se comportam como 

u(x,t) « - =Â ó v ( X 1 , quando 1 3> 1, (1-H) 

onde 

l> + 1 _____ ( /' • 1 2 



= \l-^r exp x J (L12) 

e Af = J R f(y)dy (veja o Teorema 12. 111 no Capítulo EJ) . 

Braga et al. [H] conjecturaram também que qualquer perturbação não-linear irrelevante (isto é, 
com i]f > 0) está na mesma classe de universalidade do problema linear, caracterizada pelos 
expoentes críticos a = (3 = (p + l)/2 e pela distribuição gaussiana 4> p dada em (|1.12J1 - veja, no 
Capítulo |2J a Conjectura 12. 121 

1.4 A técnica de homogeneização 

O Problema ^ foi apresentado como um modelo aproximado de um meio composto e condutor, 
caracterizado por uma estrutura microscópica periódica, ou seja, um meio não-homogêneo cujo 
coeficiente de difusão K(x) = 1 + fig{x) é periódico de período 2ir. E razoável supor que as 
variações locais de K(x) não são perceptíveis se observarmos esse meio a uma escala L ^> 2tt, 
de forma que, em grandes escalas, a heterogeneidade do meio seja cada vez menos perceptível, 
se aproximando de um meio homogéneo (isto é, com K constante). Por isso, quando nos con- 
centramos apenas no termo K(x) do Problema^ temos um problema em duas escalas. A teoria 
de homogeneização [21 1201 ES] para o Problema ^ (com A = 0, isto é, sem a parte não-linear) dá 
uma resposta rigorosa ao argumento heurístico aqui apresentado. Mais que isso, nos diz que a 
solução u se comporta como a solução de 

u(x,l) = f(x), 
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onde o = (K x ) sendo que (•) denota a média sobre o período de uma função periódica 
qualquer (essa notação será mantida em todo este texto). 

A substituição de K(x) pela sua média harmónica no caso unidimensional possui uma explicação 
heurística simples. Primeiro note que o coeficiente de difusão representa a facilidade de propa- 
gação (de calor, do local mais quente para o mais frio; de partículas, do local de maior para 
o de menor concentração; ou de cargas elétricas, do local de maior para o de menor potencial 
elétrico) enquanto sua recíproca representa dificuldade. Agora considere, por exemplo, um sis- 
tema elétrico constituído de diferentes condutores em série. Sabemos, da teoria elementar de 
circuitos elétricos [Hj, que se quisermos substituir esses condutores por outros idênticos entre si 
(isto é, com uma condutância homogénea), então o valor dessa condutância deve ser dado pela 
média harmónica das condutâncias originais. O mesmo fenómeno ocorre na difusão de calor: 
podemos imaginar o meio condutor não-homogêneo como sendo uma superposição em série de 
inúmeros condutores de espessura Ax e, se quisermos substituí-los por um condutor somente (ou 
vários com mesma espessura Ax e idênticos entre si), temos que considerar a média harmónica 
desses condutores infinitesimais. 

No Capítulo El desta dissertação vamos estudar rigorosamente alguns problemas de homogenei- 
zação. Em particular, o Corolário 15. 121 justifica a substituição de K(x) por a em (jl,13j) . 

A técnica de homogeneização foi desenvolvida para estudar o comportamento assintótico de 
soluções de equações diferenciais parciais com coeficientes periódicos. Problemas dessa natureza 
têm duas escalas de comprimento típicas: a escala do meio difusivo, associada ao período do 
coeficiente de difusão, e a escala de observação, associada ao ponto espacial e ao tempo no qual 
se fazem medidas experimentais. Por outro lado, a técnica do grupo de renormalização jHJ 
foi desenvolvida para tratar problemas com um número infinito de escalas, como por exemplo 
a equação de difusão não-linear (Problema ^ com /i = 0). Se substituirmos K(x) por o no 
Problema ^ podemos aplicar o teorema de Bricmont et al. [7] mencionado acima. Com esse 
argumento Braga et al. jlj conjecturaram que o mesmo resultado continua válido sem a hipótese 
de homogeneidade para o coeficiente de difusão (/i = 0), exceto porque a função de perfil é 
renormalizada. Em outras palavras, se F pode ser escrita na forma F(u,u x ,u xx ) = u a u b x u xx , se 
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dp > O e se A é suficientemente pequeno, então 



u(x, t) 



Vi 



A 




Vi 



x 



) 



quando t 1, 



(1.14) 



onde 



(j) a {x) 



y/Ana 



1 




) 



(1.15) 



Para uma formulação precisa de ()1.14j) . veja a Conjectura 12.91 no Capítulo El 

O termo fig(x) tem representado um obstáculo para o estudo rigoroso do Problema ^usando-se o 
grupo de renormalização analítico. Até onde sabemos, esta é a primeira vez em que se estabelece 
uma conexão entre a técnica de homogeneização e o grupo de renormalização para EDP's (no 
contexto de física estatística, veja o trabalho de Naddaf e Spencer |2~3]). 

Duro e Zuazua [TT] consideraram a equação 

u t -dw(a(x)Vu) =d-V(\u\ q - 1 u) em M 7V x(0,cx)), 

onde d G M. N e a(x) é uma matriz suave, simétrica, periódica e uniformemente elíptica. Usando a 
teoria de homogeneização eles provaram, dentre outros resultados, que a classe de universalidade 
do comportamento auto-similar da equação linear (d = 0) inclui perturbações com q > 1 + l/N. 
Apesar de não provarmos rigorosamente, nossos resultados numéricos indicam que o resultado 
de Duro e Zuazua deve continuar válido para outras perturbações não-lineares 4 . 

1.5 O grupo de renormalização numérico 

Partindo exatamente das mesmas idéias em que se baseiam o RG formal e o RG analítico, surge 
naturalmente uma versão numérica do mesmo procedimento, o chamado grupo de renormalização 
numérico, ou simplesmente RG numérico. Até onde sabemos, a primeira versão deste método 
foi proposta por Chen e Goldenfeld em 1995 [S]. O método numérico consiste, assim como no 

caso analítico, em evoluir a condição inicial segundo o fluxo dado pela EDP e depois fazer uma 

4 Chamamos atenção para uma sutil diferença entre a natureza das equações. O problema estudado por Duro 
c Zuazua, quando restrito ao caso unidimensional, inclui um termo da forma (K(x)u x )x enquanto o problema que 
estudamos computacionalmente neste trabalho tem o termo da forma (K (x)u xx ) . 
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mudança de escalas. A primeira vantagem do RG numérico, obviamente, é que não precisamos 
saber resolver a equação analiticamente ou restringir nossas hipóteses aos dados e parâmetros 
"suficientemente pequenos", uma vez que a mesma será resolvida numericamente. Outra difer- 
ença está no cálculo dos expoentes «e/3: no algoritmo numérico que vamos utilizar, os expoentes 
críticos são calculados dinamicamente, ao fim de cada iteração do algoritmo, mesmo que não se- 
jam conhecidos a priori 5 . Por último, o RG numérico permite calcular prefatores nos casos em 
que os teoremas se limitam a garantir sua existência. 

Nesta dissertação vamos apresentar a versão do RG numérico que foi implementada e utilizada 
nos trabalhos feitos em colaboração com G. Braga, F. Furtado e J. Moreira. Nesses trabalhos 
estudamos numericamente o comportamento assintótico de soluções dos problemas já apresenta- 
dos e determinamos os expoentes críticos, pré-fatores e funções de perfil, assim como tentamos 
caracterizar suas classes de universalidade. Essa versão do RG numérico, que considera a renor- 
malização das equações e permite o cálculo dos prefatores, além de uma malha discreta fixa, foi 
introduzida originalmente por G. Braga, F. Furtado e V. Isaia - veja a tese de doutorado de 
Isaia [TH] . 

Os problemas estudados nesta dissertação são todos em apenas uma dimensão espacial. Ressalta- 
mos, entretanto, que tal restrição não é, de forma alguma, uma restrição intrínseca ao método 
do grupo de renormalização numérico. Veja, novamente, a tese de doutorado de Isaia ^1 para 
o estudo de problemas com mais dimensções espaciais. 

Para o Problema Q a verificação do comportamento assintótico (|1.14jl . incluindo o estudo de 
um caso particular de perturbações relevantes, bem como a descrição e a validação do método 
numérico utilizado, foram publicados na revista Multiscale Modeling and Simulation da 
SIAM. 

Para o Problema |21 além de verificarmos o comportamento (jl.llj) . fizemos o estudo do compor- 
tamento próximo da linha que separa perturbações relevantes e irrelevantes. Neste caso, nos 

5 Isso é particularmente útil no caso de equações cujos expoentes críticos não podem ser obtidos a partir de 
simetrias da EDP, como no caso da equação de Barenblatt. Neste caso dizemos que o problema possui expoentes 
críticos anómalos. 



19 



restringimos ao caso particular da perturbação nao-linear 



XF = -u\ 



onde a G (l,oo), e trabalhamos no espaço de parâmetros (p,a). A condição para que o termo 
— u a seja marginal (r)p = 0) se reduz a 



cujo gráfico no plano (p,a) é o que chamamos de curva crítica (veja a Figura H.19|) . Estudando 
numericamente o comportamento de soluções próximas a essa curva e analisando o expoente 
crítico a, verificamos que a mesma é de fato a interface entre as regiões relevante e irrelevante e 
dizemos que o sistema sofre uma transição de fase sobre esta curva. Estes resultados estão sendo 
redigidos na forma de artigo [3] a ser submetido para publicação. 

Partes do conteúdo desta dissertação e da dissertação de mestrado de Moreira [22] estão também 
publicadas nas notas |3] do mini-curso oferecido no 57° Seminário Brasileiro de Análise. 

Todos os casos aqui mencionados são discutidos no Capítulo 121 e os resultados numéricos são 
apresentados no Capítulo |U 

1.6 Considerações finais 

Este trabalho está dividido da seguinte maneira: no Capítulo 121 estudamos formalmente os Pro- 
blemas ^e|21 através do grupo de renormalização. Antes disso, as idéias e a motivação do método 
serão introduzidas para a equação do calor. No Capítulo El descrevemos uma versão numérica do 
grupo de renormalização. O método numérico consiste na integração da equação por um curto 
intervalo de tempo, seguido de um reescalonamento. O mesmo fornece um quadro detalhado do 
regime assintótico: função perfil, expoentes críticos, pré-fatores e parâmetros efetivos. O método 
é aplicável a problemas com soluções assintoticamente auto-similares, inclusive no caso crítico 
em que aparece correção logarítmica (decaimento com (t\nt) a ). O Capítulo 0] contém, principal- 
mente, a verificação numérica das conjecturas e dos teoremas apresentados no Capítulo 121 com 
ênfase na verificação de 1)1.14)1 e (jl.llj) . além da transição de fase em (jl.lfij) . No Capítulo El vamos 



a = 



p + 3 
p+V 
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apresentar a teoria de homogeneização para equações lineares. Um dos objetivos é obter resulta- 
dos rigorosos que justifiquem a substituição do coeficiente não-homogêneo por um homogéneo no 
estudo do Problema ^ Pretendemos apresentar ao leitor um texto auto-contido sobre a teoria de 
homogeneização, para equações elípticas e parabólicas lineares, em espaços de Sobolev. O leitor 
deve estar familiarizado com a Teoria de Medida e Integração de Lebesgue e, caso não tenha feito 
um curso de Análise Funcional, deverá ter um livro introdutório à mão (por exemplo, [Hl]) para 
buscar alguns conceitos e resultados básicos. 

Listamos aqui alguns problemas que não são tratados neste texto para os quais acreditamos 
ser possível - e interessante - fazer um estudo numérico cuidadoso. Alguns deles chegaram a 
ser estudados pelo autor e colaboradores, mas não foram concluídos. No Problema ^ pode-se 
substituir o termo difusivo por um que seja conservativo, a saber, u t = {K(x)u x ) x . O estudo da 
equação Ut = [K(x)(u m ) x ] x , motivada pela difusão em meios porosos não-homogêneos, é também 
de grande interesse; o autor desconhece qualquer estudo formal, numérico ou analítico acerca da 
homogeneização desta equação. Também o estudo da transição de fase no espaço de parâmetros 
(m, a) da equação u t = {u m ) xx —u a é um representante de todo um conjunto de problemas (em que 
se enquadra o Problema 121 desta dissertação) em que um parâmetro do termo difusivo determina 
o expoente crítico de decaimento e há uma competição entre o regime linear e não-linear, dando 
origem a uma transição de fase. Em todos os problemas em que ocorre essa transição de fase, 
mesmo que com apenas um parâmetro (como u t = u xx — u a ), o estudo numérico fino das correções 
logarítmicas que surgem nos pontos críticos é também questão interessante na opinião do autor. 
Em outra direção ortogonal, podemos pensar em adaptar as idéias do método para o estudo de 
equações (não necessariamente hiperbólicas) para as quais a solução se comporta assintoticamente 
como uma "onda viajante". Para este último caso o autor chegou a implementar uma versão 
que forneceu resultados preliminares bastante satisfatórios. Por último, a formulação de um 
método mais geral, que seja apropriado para estudar problemas que apresentam decaimento, 
espalhamento e translação em seu comportamento assintótico, é questão de grande importância 
e alguns dos colaboradores do autor estão trabalhando nessa direção. 

Por último, observamos que decidimos definir 0*, CT e cf> p por (|1.9|) . (|1.15|) e í)1.12jl . respectiva- 
mente, para manter o texto o mais limpo possível. Com esta notação, o sub-índice não é um 
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parâmetro propriamente dito; por exemplo, = — (j) p se a = 1 e p = 0; assim, nao faz sentido 
escrever (f)^u. 
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Capítulo 2 

O Grupo de Renormalização Formal e 
Analítico 

Neste capítulo vamos estudar, usando o método do Grupo de Renormalização - RG, os problemas 
apresentados na Introdução. As ideias do método, assim como sua motivação, são introduzidas 
de forma natural para a equação do calor, devido às suas simetrias. Por essa razão a Seção 12.11 
é dedicada ao estudo desta equação. Nessa seção nós introduzimos a operação por mudança 
de escalas e mostramos que a equação do calor tem uma solução simétrica com respeito a esta 
operação. Essa solução é utilizada para provarmos a Proposição 12 .2\ sobre o comportamento 
assintótico de soluções do problema de valor inicial associado. 

Motivados pela operação de mudança de escalas, na Seção 12.21 nós definimos o operador RG. 
Vamos utilizar alguns resultados básicos sobre a transformada de Fourier aplicada à equação 
do calor. Nossa preocupação principal não está no preenchimento detalhado das definições e 
demonstrações, mas na idéia do RG como uma ferramenta para se provar teoremas a respeito 
do comportamento assintótico de soluções de equações difusivas. Por essa mesma razão, o leitor 
que não esteja familiarizado com esses resultados pode ler o texto sem se preocupar com tais 
detalhes. O nosso objetivo é apresentar uma outra prova da Proposição I2.2| utilizando, para isto, 
o reescalonamento da solução do pvi iterativamente. Os ingredientes básicos para a realização 
desse procedimento são: a Propriedade de Semi grupo 12.31 e o Lema da Contracão 12. 41 

Na Seção 12.31 os Problemas Q e 121 da Introdução serão considerados sob o ponto de vista formal 
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e heurístico. Observamos que a aplicação do método descrito acima depende fortemente da 
invariância da equação do calor por uma mudança de escalas. Contudo, em geral as equações 
não são invariantes sob este reescalonamento. Nesse caso, para analisar os Problemas Q e 121 sob o 
ponto de vista do RG, nós consideraremos uma equação diferente em cada escala. Dizemos que 
a equação foi renormalizada. No processo de renormalização das equações, esperamos que fique 
latente a simplicidade e o poder da mudança de escalas: é possível determinar toda uma classe 
de perturbações não-lineares F, para as quais o comportamento assintótico das soluções continua 
sendo o mesmo da equação linear. 

2.1 A equação do calor e mudança de escalas 

Nesta seção vamos considerar o problema de valor inicial para a equação do calor: 



onde Co(R) é o conjunto das funções contínuas de suporte compacto em R. 

Enunciamos, a seguir, o teorema que nos garante a existência (global) e unicidade do pvi (J2.1J) . 
como também nos fornece uma representação integral para a solução. Esse teorema pode ser for- 
mulado sob condições mais fracas, mas para os nossos propósitos é suficiente o caso em que o dado 
inicial seja contínuo de suporte compacto. A prova da proposição abaixo pode ser encontrada 
em UB]. 

Proposição 2.1 problema de valor inicial Wl\ com dado inicial f G Co(R), possui uma única 
solução u(-,t) G C°°(R), t > 0. Essa solução é dada pela fórmula 



Ut = u xx , x G R, t > 

U(;0) = f, /GCo(R), 






para qualquer t > e u(x,t) — * f(x) quando t J. 0. 



Seja L um número real estritamente positivo e façamos a seguinte mudança de escalas: 

x i— ► Lx t i— ► L 2 t u i— > Lu. 
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(2.3) 



Chamemos de v(x,t) à função u(x,t) devidamente reescalonada, isto é, 

v(x,t) = Lu(Lx,L 2 t). (2.4) 

Temos que v t = L 3 u s e v xx = L 3 u yy , onde s = L 2 t e y = Lx, ou seja, u t = u xx se e somente 
se v t = v xx . Dizemos que a equação do calor é invariante pela mudança de escalas (|2.4|) . para 
qualquer L > 0. 

Vamos procurar soluções que satisfaçam u(x,t) = Lu(Lx,L 2 t). Substituindo formalmente L por 
l/Vi, temos u(x,t) = ^ u(x/Vt, 1), o que nos motiva a procurar soluções da forma 

" ( *< í) = Mti)' (2 ' 5) 

Qualquer função n que se escreve na forma (J2.5|) é obviamente invariante pela mudança de 
escalas ()2.4|) . Para determinar a função </>(■), nós substituímos ()2.5|) em í)2.1|) . obtendo uma 
equação diferencial ordinária para </>(£) : 

= 0" + 1^' + Í0 = (V + . (2.6) 



2" r 2 r V 2, 

A equação acima admite soluções da forma C exp(— ^ 2 /4) e, escolhendo a constante C de modo 
que L(j>(Ç)d£ = 1, obtemos a distribuição gaussiana 0* dada por (|1.9|) . Substituindo por 0* 
em f)2.5j) . temos uma solução invariante pela mudança de escalas, que satisfaz a seguinte relação: 
yjiu(y/i x,t) = 4>*{x). A proposição abaixo mostra como essa relação está relacionada com o 
comportamento assintótico da solução geral. 

Proposição 2.2 Seja u(x,t) a solução da equação do calor com dado inicial f G Cq(M). Então 

Vi u(Vi x, t) : — > A f( f)*(x), V x G IR, (2.7) 

onde A f = J R f(y)dy. 



Prova: Substituindo x por x\ft e multiplicando (|2.2|) por Vi, obtemos 

/— /— 1 / —(zVt-y) 2 i 

Vt u(Vt x,t) - Af<f)*(x) = —== / e b f(y)dy-(f)*(x) / f(y)dy 



4tt 
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= J(^e^^-Mx)jf(y)dy (2.8) 

= Jj,(x)íe ? ^ Z -l ) jf(y)dy. (2.9) 

Primeiro note que o integrando de ()2.8|) é dominado em módulo por \f(y)\/y/n, que é integrável. 
Note também que o integrando de ()2.9|) converge a zero quando t — > oo, isso para todo y G 
M (podemos considerar x como um parâmetro fixo). Aplica-se o Teorema da Convergência 
Dominada de Lebesgue para se obter (|2.7|) . | — | 



2.2 A transformação do grupo de renormalização 

O limite na Proposição 12 . 21 diz que soluções da equação do calor têm o comportamento assintótico 
ditado pela solução invariante por mudança de escalas. Vimos anteriormente que a equação 
do calor é invariante por essa mudança. Isso nos motiva a definir o operador do grupo de 
renormalização (veja jElIZl)- Esse operador age num espaço de funções, mais especificamente no 
espaço dos dados iniciais, e o mesmo consiste em evoluir, no tempo, o dado inicial f(x), seguindo 
o fluxo dado pela equação do calor (essa evolução é realizada do tempo t — 1 ao tempo t = L 2 ) 
e, em seguida, mudamos a escala do domínio e do contra-domínio do dado evoluído. Denotando 
esse operador por R L , onde L > 1 é o fator de escala 1 , podemos descrever essa sequência de 
operações pelo diagrama: 

/ = U f {; 1) ' U f {; L 2 ) h-> LUf(L ; L 2 ) = R L f, (2.10) 

onde Uf denota a solução da equação do calor com condição inicial u(x, 1) = f(x). De outra 
forma: 

(R L f)(x) = Lu f (Lx,L 2 ). (2.11) 

E claro que Rl é um operador linear, uma vez que a equação do calor é linear e homogénea. 

Vamos estudar como Rl atua conjugado com a transformada de Fourier (veja a definição e 

principais propriedades dessa transformada no Apêndice 16. lj) . 
1 Daqui cm diante vamos sempre assumir que L > 1. 
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No caso da equação do calor, para a qual já sabemos da existência e unicidade de uma solução 
global, a Propriedade de Semigrupo, que enunciamos e provamos abaixo, segue trivialmente da 
definição de R^. Entretanto, vamos usar a transformada de Fourier para ilustrar um outro 
aspecto do RG analítico: podemos provar a existência da solução global do problema a partir da 
existência local. 

Primeiramente, se u denota a única solução da equação do calor com condição inicial u(-, 0) = /, 
e Uf denota a única solução com condição inicial «/(■, 1) = /, então temos 

ú(k,t) = e~ kH f(k) e ú f (k,t) = e- k2{t - l) f(k). (2.12) 



Para uma prova deste fato, veja [TH] - Aplicando ()6.2|) a (J2TTJ), e substituindo (|2.12J) . temos 

(S7)(fc) = (Lu^L^k) = Ú f (k/L,L 2 ) = e-W-Vf(k/L). (2.13) 
Uma propriedade fundamental do operador R L é a seguinte: 

Teorema 2.3 (Propriedade de Semigrupo) Sejam L\ > 1 e L 2 > 1. Então 

R Ll oR L2 = R LlL2 . (2.14) 

Prova: Basta desenvolver (|2.13j) para mostrar que (Rl 1 Rl 2 Í) — (Rl 1 l 2 Í) e ° resultado segue, 
como vemos: 

(RZR7J) (*) = e-W^M-QíSZjHk/Ia) 

= e~ k2{1 ~^ ) e~ k2( ^^ ) f^/LM 

= (íCj) (*). 

A prova está concluída. ^ 

Também podemos "diagonalizar" o operador Rl, que tem os polinómios de Hermite com peso 
gaussiano como auto-vetores no espaço de Fourier. 

A transformada de Fourier da gaussiana 0* (x) dada por (JTJJ) é fa{k) = e~ k . De temos 

(Rrf.)(k) = e- k ^-^e-Ti = fak), (2.15) 
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ou seja, a distribuição gaussiana 0* é um ponto fixo do RG. Da mesma forma, se f n (x) é a 
função cuja transformada de Fourier é f n (k) = k n e~ k2 (esses são, essencialmente, os polinómios 
de Hermite com peso gaussiano), então f n é um autovetor de Rl com autovalor (1/ ' L) n . Dizemos 
que a gaussiana <fi*(x) está na direção marginal e que os outros autovetores estão na direção 
irrelevante. O subespaço gerado pelos autovetores irrelevantes também é chamado de variedade 
estável. 



Definimos o espaço £> 4 por 



B A = [f(x) G L 2 (R) : / G C\ 
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< + OO 



onde 



sup 

fcgR 



;i + ^ 4 )(i/(£0i + i/»i) 



Pode-se mostrar que B4 é um espaço de Banach, que os elementos de B4 são funções contínuas 
e que a norma || • ||,g 4 domina as normas L X (R), L 2 (M) e L°°(IR) (veja, por exemplo, a tese de 
mestrado de Moreira [22] ou as notas de mini-curso 0). 

Teorema 2.4 (Lema de Contração) 5e / G B4 e u(x,t) denota a solução da equação do calor 
tal que w(-, 0) = /, então w(-, t) G £>4 V t > 0. Em particular, RlB^ C B4. 



Além disso, se /(O) = 0, vale 



c 

~L 



S 4 ) 



onde C = C(L) > é decrescente em relação a L e não depende de f . 



Prova: Começaremos pela segunda parte. A equação do calor pode ser vista como um caso partic- 
ular do Problema El com p = X = e K(t) = 1. A prova que apresentamos aqui é uma adaptação 
para este caso particular de um prova mais geral - veja o Lema de Contração em jH]. A partir de 
(j2H temos RTf(k) = e "(r) 2 ( i2 -i) f(k/ X) e R^~f(k) = -(L 2 - \){2k/ L 2 )e-^) 2 ^ 2 ^) f{k/ L) + 
L^ e ~^) 2 (L 2 -i) f{k/L). Logo, vale que 



\R L f(k)\ + \R L f (jOj 

(*) 2(i 2_l) 



< 



f(k/L) + 2kf(k/L) 



-f\k/L) 
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k/L 

/(O) + / f\s)ds 



o 



/ lfc/i ' fx^l 

O 1 T S 



Valem as estimativas f'(k/L) < ||/||b 4 e 
l/>/i)l 

uma vez que /(O) = 0. Assim, temos 

(1 + + |^/ (*)|) < (1 + fe)(1 + |A:|+ / |A: 1)6 ||/|k. 

Definindo C como 

C(L) = sup {(1 + k l ) (1 + \k\ + 2|£; 2 |) e-tr) 2 ^ 2 " 1 )) < +oo, L > 1, 

temos que ||-Rl/||b 4 < ^||/||b 4 - Além disso, C(L) é positivo e decrescente, logo C/L é estrita- 
mente decrescente e lim^oo C/L = 0. 

A primeira parte pode ser mostrada de forma idêntica, exceto porque ()2.16j) não vale em geral, o 
que não é problema, pois a exponencial de (j2.12|) domina qualquer polinómio e com isso podemos 
mostrar que u(-, t) G £>4 V t > 0. 

Observação: Segue do Teorema 12 .41 que existe L > 1 tal que Rl é uma contração em {/ G £>4 : 
/(O) = 0} para qualquer L > L . 

Corolário 2.5 (Convergência) Se L > L e / G i3 4 entôo 

iÇ/ ^ Af0*> (2.17) 

n— >oo 

onde A/ = J R f(y)dy. 

Prova: Defina f n = R L f e g n = f n — Afif)*, n = 0, 1,2, . . .. Pela linearidade de Rl e como 
0* = Rl<P*, temos que g n+l = f n+1 - A f (f>* = R L f n - R L A f 4>* = R L g n - É suficiente ver 
que g^(0) = \/n G N, pois, pelo Teorema 12.41 teremos g n — -> 0. Mas gò(0) = L g (y)dy = 
IrU '- A f4>*)(y)dy = f u f(y)dy-A f f u (p*(y)dy = 0. Seja n>0e suponha que gZ(0) = 0, temos, 
por (|2.13|) . que <^+i(0) = RLg n {fy = 0, o que, por indução, conclui esta prova. |— I 

O corolário abaixo traz um resultado similar à Proposição 12.21 
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Corolário 2.6 Seja u(x,t) a solução da equação do calor com dado inicial u{- , 0) = / G £>4. 
Então vale 

^ft~ n u[y[T n x, t n ) A f <t>*(x) (2.18) 

n-^oo 

para qualquer sequência t n = L 2n com L > Lq, onde Af = j R f(y)dy. Em particular, temos 
convergência uniforme e convergência nas topologias de L 1 (M) e L 2 (R). 

Prova: Como a equação que estamos considerando tem sua condição inicial no tempo í = 0eo 
grupo de renormalização foi definido para soluções com condição inicial em t — 1, definimos a 
função g como g = u(-, 1). Segue de (|2.12|) que g(0) = /(O) e, portanto, J R g(y)dy = f R f(y)dy = 
A f . 

Agora, veja que da Propriedade de Semigrupo obtém-se i?£ = Rl« e, pela definição de Rl (|2.11J) . 
temos 

R n L g{x) = R L ng(x) = L n u{L n x,L 2n ) 

Segue do Corolário 12.51 que 

L n u(L n x,L 2n ) AMx), 

n— >oo 

concluindo a demonstração. |— | 

Observação: O limite pode ser mostrado para t — > oo em geral, preenchendo- se os intervalos 
[L 2n , L 2n+2 }. Entretanto, esse último passo foge ao objetivo desta dissertação, pois é uma parte 
mais técnica e que ilustra pouco da utilidade do RG. Além disso, estamos considerando sistemas 
cujas soluções são "bem comportadas" e é pouco razoável a possibilidade de tal limite valer para 
toda subseqúência da forma acima mas não valer para t — > oo em geral. Por último, no algoritmo 
numérico que vamos utilizar neste trabalho, somente essas subseqúências serão consideradas. O 
leitor interessado pode consultar |3J 1221 • 

2.3 Renormalização das equações 

Na seção anterior definimos o grupo de renormalização motivados pelo fato de a equação do calor 
ser invariante pela mudança de escalas (J2.4)) . O que está por trás do fato de o método determinar 
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o comportamento assintótico do pvi é a Propriedade de Semigrupo (para a equação do calor, veja 
como a demonstração do Corolário 12.61 utiliza o Teorema 12 .'d\i . Na verdade, a transformação do 
RG aplicada à equação do calor é um caso particular da aplicação desse método. No caso mais 
geral (veja [SJ Ej), a equação dada não é necessariamente invariante pela ação das operações (|2.3|) . 
Por isso temos que ampliar o conceito do RG, fazendo o mesmo agir também sobre a equação 
diferencial dada. O que se faz é definir a solução reescalonada v como em í)2.4j] e considerar a 
equação satisfeita por v. Assim, o esquema para o caso geral é: 

eq eq 1 eq 2 eq n , nin v 

/o /l J2 Jn 

onde o n-ésimo mapeamento evolui o n-ésimo dado inicial /„ pela n-ésima equação eq n , com 
n = 0,1,2,.... Portanto a transformação do grupo de renormalização terá, no caso geral, duas 
componentes: uma atuando sobre o espaço dos dados iniciais e a outra atuando sobre a equação 
diferencial em questão. Com isso mantemos a propriedade de semigrupo, fato essencial para que 
possamos fazer a evolução temporal do dado inicial usando o RG. 

Seja u(x,t), iel,í>l, solução global 2 de um dos dois problemas apresentados no Capítulo[T] 
Ao invés de tomarmos o limite (|l.fi|) de uma só vez, o mesmo será fruto de um processo iterativo. 
Para um L > 1 fixo e expoentes a > e (3 > 0, definimos 3 uma sequência (tt n )neN de funções, 
por uq =W6, para n = 1, 2, 3, . . . 

u n {x, t) = L 2a u n ^(L 2l3 x, L 2 t). (2.20) 

Diremos que (u n ) é uma sequência de soluções reescalonadas. A partir da definição ()2.20|) e 
sabendo-se qual equação diferencial é satisfeita por u n _i, é um exercício simples de cálculo obter 
a equação diferencial satisfeita por u n . A estas equações chamaremos sequência de equações 
renormalizadas. 

Vamos definir f n por 

f n (x) = u n (x,l). (2.21) 

2 Não estamos preocupados neste momento com a questão da existência global de tal solução, mas simplesmente 
assumindo este fato. 

3 A forma mais "correta" seria definir u n (x,t) — L a u n ^i{L^x,Lt). A substituição de L por L 2 neste capítulo 
pode ser vista como uma questão simplesmente estética. O RG numérico, que veremos no Capítulo utiliza L 
como fator de escala ao invés de L 2 . 
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Segue de QPIJ) e (gHH) que 

f n (x)=L 2a u n _ 1 (L 2 ?x,L 2 ). (2.22) 

Observe que / n+ i pode ser obtida evoluindo-se f n de í = 1 até í = L 2 , de acordo com o fluxo 
que determina u n (-,t), e aplicando-se as mudanças de escalas x i— »• L 2 @x, u i— » L 2a u. Observe 
também que uma consequência imediata de (|2.20|) é 

f n (x) = u n (x, 1) = L 2na u (L 2nf3 x, L 2n ) . (2.23) 

Ressaltamos a importância da relação acima. Suponha que, para o par de expoentes a e (3 
dados, sejamos capazes de mostrar que f n — > <p, onde (f) : R — > R não é identicamente nula. Neste 
caso, teremos provado o limite (|1.6|) para a subseqúência t n = L 2n . Para provar o mesmo limite 
para qualquer sequencia t — > oo de forma rigorosa, ficaria faltando então preencher os intervalos 
[L 2n , L 2n+2 } - veja a observação após o Corolário 12.61 Essa é, essencialmente, a forma com que 
a técnica do RG analítico é aplicada. A seguir, aplicaremos as idéias explicadas neste parágrafo 
para o problema de difusão não-linear. Os Problemas Q e El serão considerados nas subseções 
seguintes. 

Consideremos o problema de difusão não-linear (Problema com \x = 0). Vamos assumir que 
a perturbação é da forma F = u a u b x u c xx . Se u é solução global do problema, então a n-ésima 
equação renormalizada será aquela satisfeita por u n definida em (|2.20|) . dada por 

u t = XnU xx + X n u a u b x u c xx , (2.24) 

onde 

Xn = (£1-20)2" e \ n = ( L l-(6+2c)/3-(a+6 + c-l)a )2 n A _ 

Se estamos interessados em caracterizar os problemas cujas soluções se comportem como no caso 
da equação do calor, esperamos algum tipo de convergência e por isso não é interessante que o 
termo \n tenda a ou oo quando n cresce. Para isso, devemos ter (3=1/2. Como a equação do 
calor é conservativa 4 , devemos ter a = (3. Substituindo a = 1/2 — (3 temos que 

X„ = l e \ n = (L~ dF ) n A, 

4 Dizer que a equação é conservativa significa que suas soluções satisfazem J R u(x,t)dx = J m f(y)aly V t > 1 
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de forma que a parte nao-linear vai tender a se anular, permanecer constante, ou tender a infinito, 
se dp = a + 26 + 3c — 3 for positivo, nulo ou negativo, respectivamente. 

Ainda que não constitua uma prova rigorosa, o argumento acima permite classificar, mediante 
um exercício simples de derivação, uma gama de perturbações não-lineares da forma F = u a u b x u c xx 
que estão na mesma classe de universalidade 5 da equação linear, caracterizada pelos expoentes 
críticos a e j3. 

Com essa abordagem, Bricmont et al. jT] provaram o teorema abaixo 6 . 

Teorema 2.7 (Bricmont et al., 1994) Considere o Problema^ com \x = e condição inicial 
u(-,0) = f G B4. Suponha que a perturbação F possa ser escrita na forma F(u,u x ,u xx ) = 
u a u b x u c xx . Nessas condições, se a + 26 + 3c — 3 > e A é suficientemente pequeno, então existe 
AeR tal que 

Vi u(y/i x,t) A<fi*(x), 

n— >oo 

onde u é a solução do problema de valor inicial e 0* é definido em M.ty) . 

Em outras palavras, Bricmont et al. mostraram que u é assintoticamente auto-similar da forma 

Definição 2.8 (A transformação do grupo de renormalização) Paran = 0, 1,2, . . defin- 
imos o operador R L n como 

R L ,nf = Lu f (L-,L 2 ), 

onde Uf(x,t) denota a solução da n-ésima equação renormalizada \2. 24\ ) com condição inicial 
u f {x,l) = f{x). 

De acordo com a Definição 12.81 e de acordo com as equações (|2.2U|) , (|2.21|) e (|2.22|) , temos que 

fn+l(x) = Lu n {Lx,L 2 ) = R L ,nfn- 
5 Veja a discussão sobre classes de universalidade na Introdução, página [TUI 

6 Um material acessível e minuncioso em que estudam-se perturbações da forma F = u a , a = 4, 5, 6, . . . , é a tese 
de mestrado de Moreira [221 ■ O resultado de Bricmont et al. é mais geral do que o Teorema 12 . 71 que enunciamos 
aqui. 
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Usando a identidade acima e denotando por R\ a composição Rl,u-i ° Rl,u-2 o • • • o R L1 o Rl,o, 
segue de (I2.23j) que 

f n {x) = L n u(L n x, L 2n ) = R L n fi f(x) = R n J(x), n = 1, 2, . . . , (2.25) 
ao que chamaremos de "propriedade de semigrupo". 

2.3.1 Difusão com coeficiente periódico 

Considere o Problema [T] com F = u a u b x u c xx . Seguindo o mesmo procedimento feito para a equação 
de difusão não-linear acima, temos que a n-ésima equação renormalizada será 

Ut = Xn [1 + iig(u) n x)]u xx + A n u a u b x u c XXJ (2.26) 

onde 



n = (L 1_2 ^) 2n , Wn = (L /3 ) 2n , X n = A (L 1 -( b+2c ^( a+b+c - 1 ) a ) 2 ". 



X 



Observe que estamos no mesmo caso da difusão não-linear exceto pela adição do termo fig(x). 
No caso de /i = 7^ A, o Teorema 12.71 pode ser aplicado para garantir que 

u(x,t) w -A ^* ( 4=) ■ ( 2 - 27 ) 



Gostaríamos de entender como o termo ng(x) afeta a classe de universalidade do comportamento 
assintótico auto-similar dado por (|2.27[) . Para isso, é natural escolher novamente a — [3 — 1/2. 
Nesse caso temos 

X„ = l, co n = L n , \ n = \(L- d n n - 

Consideremos, inicialmente, o Problema 1 supondo que \x ^ = A (equação linear com coeficiente 
de difusão periódico). Repare que, a cada renormalização, temos uma equação similar, porém 
com o termo fig(x) oscilando cada vez mais rapidamente. A teoria de homogeneização, que 
estudaremos no Capítulo |3J nos diz que, para grandes escalas, as soluções desse problema - 
dito não-homogêneo, pois reflete as propriedades de um meio não-homogêneo - se aproximam 
da solução de uma equação com coeficiente constante, chamada equação homogeneizada. Nesse 
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caso, conforme a discussão feita no Capítulo ^(e que será vista com mais detalhes no Capítulo [HJ), 
este coeficiente constante é dado pela média harmónica do coeficiente 1 + [ig(x). Tome e" 1 = L n . 
Observando que e n — > 0, vemos que o estudo do comportamento assintótico da solução desse 

n— >oo 

problema é equivalente ao estudo, quando e — > 0, da solução da equação 

u t =[l + jigie^x^u^. 

De acordo com a teoria de homogeneização, a solução da equação acima se aproxima, quando 
e — > 0, de um múltiplo da solução de 

ut = o-u xx , (2.28) 
onde cr = ((1 + fig)^ 1 )' 1 é a média harmónica do coeficiente de difusão. 

Consideremos, agora, o caso geral para o Problema 1, isto é, suponha agora que À / ^ /i 
(equação não-linear com coeficiente de difusão periódico). A equação renormalizada é dada 
por í)2.26jl . Quando n — > oo, o termo não-linear X n F é supostamente levado a zero se d>F > 0, 
enquanto que o coeficiente periódico 1 + fig(u n x) oscila cada vez mais rapidamente, produzindo 
um resultado cada vez mais parecido com o de um coeficiente constante o = ((1 + fig)^ 1 )^ 1 . 

Somos então levados a supor que as soluções do Problema 1, dentro de uma classe de condições 
iniciais adequada e com F satisfazendo dp > 0, se comportam assintoticamente como as soluções 
de f)2.28j) . Fazendo a mudança de variáveis ihi/ y/ã, obtemos a equação do calor. Concluímos 
então que as soluções do Problema 1 se comportam da mesma forma que as soluções da equação 
de difusão não-linear, exceto pelo fator \fõ no argumento da função de perfil 0*. Com esses 
argumentos, Braga et al. [3] estabeleceram e verificaram numericamente a conjectura abaixo. 

Conjectura 2.9 (Braga et al., 2003) Considere o Problema^ Suponha que a perturbação F 
possa ser escrita na forma F(u,u x ,u xx ) = u a u b x u c xx . Se a condição inicial é dada por u(-,0) = 
f G C^°(1R), a + 26 + 3c — 3>0 e X é suficientemente pequeno, então existe uma única solução 
u(x,t), iGl,i>0. Além disso, existe AeR, tal que 

Vtu(Vix,t) — > A(j) a (x) quando t — > oo, 

onde <f) a é definido em hl.lfy) . o limite valendo em L X (IR) e L°°(R). 
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Em outras palavras, temos que u é assintoticamente auto-similar, da forma 



u(x, t) 



Vi 



A 




Vi 



x 



) 



com expoentes críticos a = (3 = \ e função de perfil (f) a . 
Um caso de perturbação relevante 

Para um caso particular de perturbações não-lineares relevantes, é possível ainda fazer o estudo 
do comportamento assintótico das soluções a partir da renormalização da equação. 

Considere o Problema ^ com XF = —u a , a G (l,+oo). O mesmo argumento que justifica a 
Conjectura 12.91 permite classificar as perturbações da forma acima como irrelevantes se a > 3. 
Pensemos agora no caso a < 3 e fixemos nossa atenção na parte não-linear do problema. A 
renormalização com a = 1/2 faz esse termo crescer. Podemos pensar que isso se deve ao fato 
de que o termo não-linear produz um decaimento mais rápido do que í" 1 / 2 - as soluções de 
u t = —u a são da forma u = A(p(x)t~ 1 ^ a ~ 1,) com qualquer <\>. Por outro lado, a parte linear não 
pode ser descartada pois ela é responsável pela difusão. Nesse caso, a renormalização correta 
é dada pelos expoentes (3 = 1/2, a = l/(a— 1). Essa escolha de (3 garante que a parte linear 
seja invariante pela renormalização. A escolha de a garante o mesmo para a parte não-linear. 
Essa nova escolha de expoentes determina um novo operador e uma nova renormalização da 
equação. Em relação a essa nova renormalização, tanto u xx quanto u a são termos marginais. 

Para o caso homogéneo (/x = 0), pode-se seguir o mesmo raciocínio que nos levou à equação 
()2.6|) para obter a equação a ser satisfeita pelo ponto fixo do novo RG. Essa equação possui uma 
solução (p a positiva, de quadrado integrável, que é única a menos de uma constante multiplicativa 
- veja j7] e referências lá citadas. Bricmont et ai. [7] provaram que se a condição inicial u(-, 0) = / 
satisfaz / G £> 4 e f(x) > V x G M, então 



Para o caso em que \x ^ 0, Braga et al. [1] verificaram numericamente a conjectura abaixo. 
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Conjectura 2.10 (Braga et al., 2003) Considere o Problema^ com condição inicial u(- , 0) = 
/ G C^°(R), / > 0. Suponha que perturbação F possa ser escrita como F(u) = u a para u > 
e que A < 0. Se a 6 (1,3) e A é suficientemente pequeno, então existe uma única solução 
u(x,t), x G R, t > 1. Além disso, 

t^u(Vix,t) — > f a ( -7=) quando t — > 00, (2.29) 

onde o - = ((1 + fig(x))' 1 )^ 1 e f a é uma função tal que 

é solução do problema com ji = 0. 



2.3.2 Difusão com coeficiente dependente do tempo 



Considere o Problema [T] com F = u a u b x u c xx . Reescrevemos a equação como: 



u t = [t p + o(t p )] u xx + A„ u a u b x u c xx . 



Neste caso, a n-ésima equação renormalizada é dada por 

o(L np t p ) 



ut = Xr 



t" 



Uxx ~t~ A n u u x u Xi 



onde 



X, 



{b+2c)0- (a+b+c- l)a \ 2n 



(2.30) 



Fazendo (3 = a = p ^-, temos 



Xn = 1, A n = A(L 



-(p+l)í? F \n 



onde 



?7i? = a + 2b + 3c 



p + 3 
p + 1' 



(2.31) 



Observe que lim^oo o(s)/s = e, portanto, 



o(L n np) 

sup v — ; — ► 0, 



Í6[1,L 2 ] 
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de forma que, para grandes valores de n, este termo se torna desprezível na equação ()2.30|) . Por 
outro lado, se r)p > 0, temos que \ n decresce rapidamente a quanto n — > oo. Assim, somos 
levados a supor que o regime assintótico das soluções do Problema 121 com r/p > 0, e cuja condição 
inicial u(-,t) — f é suave e decai suficientemente rápido, está na mesma classe de universalidade 
da solução auto-similar de 

u t = t p u xx . (2.32) 

De fato, 



u(x, t) = 




onde a = (3 = e 4> p é dada por (|1.12|) é uma solução auto-similar de (|2.32l) . 
Para o caso linear /i = 0, Braga et. al jHJ obtiveram o seguinte resultado. 

Teorema 2.11 (Braga et al., 2004) Considere o Problema^ linear, isto é, com X = 0. Se a 

condição inicial u(-,0) = / G £>4, então a solução u(x,t) satisfaz 

Vp^u(Vt^-, t) A f P (-), 

t — >oo 

onde 4> p {x) é dada por hl.lty) e Af = Lf(y)dy. 

Em outras palavras, u é assintoticamente auto-similar com expoentes a = (3 = e o prefator 
Af depende apenas da integral do dado inicial /. 

Para o caso não-linear, Braga et al. propuseram a seguinte conjectura. 

Conjectura 2.12 (Braga et al., 2004) Considere o Problema^ com F = u a u b x u c xx e defina t)f 
por \2.31)) . Então, se r/p > e A é suficientemente pequeno, continua valendo a conclusão do 
Teorema \2.11i substituindo- se Af por uma constante A, que não depende apenas de f . 

Transição de fase 

Considere a equação 

u t = t p u xx - u a , (2.33) 
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com espaço de parâmetros (p, a). 

Observe que a condição para que o termo — u a seja marginal, isto é, para que r\p dada por (j2.31j) 
se anule, se reduz a 

a = a c (p) = £±f (2.34) 
p+l 

Vamos analisar como o par de parâmetros (p, a) afeta o expoente crítico a. 

Para a > a c (p) temos que o termo — u a é supostamente irrelevante pois tende a se anular sob a 
renormalização dada pelos expoentes a = j3 = (p+l)/2. Nesse caso, o comportamento assintótico 
das soluções está na mesma classe de universalidade da equação linear e dizemos que estamos no 
regime linear. 

Para a < a c (p), temos que o termo — u a tende a crescer com as sucessivas renormalizações se os 
expoentes forem escolhidos como no parágrafo acima. Nesse caso, a escolha correta dos expoentes 
é a = l/(a — l)e(3 = (p+l)/2. O espalhamento com é causado pelo termo t p u xx e o decaimento 
com t a é determinado por — u a . O termo não- linear, que é considerado relevante, passa a ser 
marginal em relação a esses novos expoentes críticos. Nesse caso, dizemos que estamos no regime 
não-linear. 

O caso a = a c (p) é o mais delicado. Ambos os termos linear e não-linear produzem o mesmo ex- 
poente crítico de decaimento a = = ^ry- Nesse caso, dizemos que estamos no regime critico. 
No regime crítico aparece uma correção logarítmica ao decaimento assintótico das soluções. 

Podemos pensar que o termo determinante do expoente a será aquele que, sozinho, causaria o 
maior decaimento. Nesse caso, temos 

a(p, a) = max 

Observe que a superfície dada pelo gráfico de a(p, a) tem suas singularidades justamente sobre a 
curva a = a c (p). A Figura li. 191 mostra o gráfico de a c (p) separando as regiões dos regimes linear 
e não-linear. A Figura E.20I mostra o gráfico de a(p, a) com singularidades em a = a c (p). 

Conjectura 2.13 (Braga et al., 2004) Considere a equação \2. 33^ com p e [0, oo) e a E 



P + l 1 1 
\ 2 'a-lj 
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(l,oo) ; com condição inicial u(-,0) = f satisfazendo f G C^°(R) e f(x) > O V x G IR. Então 
existe uma única solução u(x,t), x G R, t > 0. 

Se a > a c {p), então existe A > O, giíe depende de p, a e f , tal que 

Vt^ U (Vt^-, t) A* &,(■), 

í^oo 

onde P e dado por M.lfy) . 

Se a = a c (p), então existe A > 0, gne depende apenas de p, tal que 

y/(t \nt)^u{V¥^-, t) 

í— >oo 

Se a < a c (p), então existe A > 0, que depende apenas de a e p, e 4> Pt a G L 1 (R) estritamente 
positiva tal que 

t^u(V^;t)-^A f (f> Pía (-). 

t— >oo 

No Capítulo El vamos fazer um estudo em torno do ponto (p, a) = (1/2,7/3), que está sobre a 
curva crítica 1)2.34)1 . Verificamos, entre outras coisas, que este é de fato um ponto de transição 
de acordo com a Conjectura 12.131 acima. Em particular, vamos variar p e a separadamente e 
calcular o expoente a numericamente, verificando a fórmula ()2.35j) . 
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Capítulo 3 

O Grupo de Renormalização Numérico 



Neste capítulo vamos descrever uma versão numérica do grupo de renormalização. O algoritmo, 
baseado nas ideias do RG analítico, é dotado de elegância e simplicidade que merecem destaque. 
Ademais, o método fornece um quadro detalhado do regime assintótico: função perfil, expoentes 
críticos, prefatores e parâmetros efetivos, sendo aplicável não só a problemas com soluções auto- 
similares (que decaem como ~ t~ a ), mas também ao caso crítico (ou marginal), em que aparecem 
correções logarítmicas ao decaimento t~ a (decaimento como ~ (ílní)~ a ). 

Por tudo isso, o RG numérico parece ser a ferramenta ideal para estudar quais aspectos do com- 
portamento assintótico são universais e quais dependem de parâmetros da EDP ou da condição 
inicial. O método também dá suporte a trabalhos rigorosos, pois os resultados obtidos podem 
encorajar tentativas de se provar teoremas ou, eventualmente, apontar para a necessidade de 
mudanças nas hipóteses. 

Pelo que sabemos, a primeira versão numérica do RG para EDP's foi proposta por Chen e 
Goldenfeld |Sj. Uma versão modificada foi introduzida por Braga, Furtado e Isaia, sendo bem 
mais rápida e eficiente que a primeira. Uma descrição detalhada desta versão pode ser encontrada 
na tese de doutorado de Isaia ^Hj- Infelizmente, parte do ganho computacional desta última não 
se aplica ao Problema ^ como veremos mais adiante. 
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O método aqui apresentado aplica-se a equações da forma 



u t = G(u,u x ,u xx ,x,t) 
u(x, 1) = f(x) 



(3.1) 



Observamos que os problemas apresentados no capítulo introdutório e que serão estudados nu- 
mericamente nesta dissertação podem ser escritos na forma acima. 

A ideia de uma sequência de mudanças de escalas e de renormalização de equações será gene- 
ralizada para sequências em que os expoentes de escala a e j3 podem ser diferentes em cada 
passo. Essa generalização será vista na Seção I3~T| juntamente outras noções preliminares de que 
vamos precisar para explicar o funcionamento do método. Cada iteração do método consiste 
basicamente nos seguintes passos: evoluímos uma dada função por um curto intervalo de tempo, 
segundo o fluxo dado pela equação (j3.1j) : mudamos as escalas do domínio e do contra-domínio 
da função obtida; renormalizamos a equação (J3.1|) . Esse algoritmo numérico será descrito na 
Seção 13.21 Outras decisões de implementação, como a discretização da equação, a escolha de 
uma escala Lea mudança de escalas no domínio da função serão discutidas na Seção 13.31 A 
validação do método numérico será feita no Capítulo |U 

3.1 Preliminares 

Seja Uq uma função real de (x,t) G IR x IR + . Para um L > 1 fixo e sequências de expoentes 
positivos (a n ) n£ N e (/3 n ) n6 N, defina a sequência de funções (u n ) ng N, indutivamente, pela seguinte 
mudança de escalas: 



Daqui em diante vamos supor que a função original uq é solução global da Eq. í)3.1j) com f = fo e 
G = Gq, de forma que (u n ) seja uma sequência de soluções reescalonadas. Neste caso, u n satisfaz 




(3.2) 



Repare que a definição ()3.2|) acima implica em 




(3.3) 
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ao seguinte problema de valor inicial renormalizado: 

i^t G n (u, u x , u xx , x, t) , 
u{x, 1) = f n {x) ' [ÓA) 

onde G n e /„ são obtidos recursivamente, como explicamos a seguir. 

Primeiramente, o operador G n pode ser obtido a partir de aplicando-se substituição de 

variáveis e a regra da cadeia 

G n (u, u x , u xx , x, t) = L an+1 G n -i (L an u,L an ^ n u x ,L an 2 n u xx , L^ n x, Li). (3-5) 

Já a condição inicial /„ pode ser obtida a partir de ] n -\ integrando-se a equação para u n _i 
no intervalo de tempo [1,L] e mudando-se as escalas do espaço e da própria função, pois, pela 
definição 1)3 .2|) . temos 

f n (x) = L a "u n - 1 {L^x,L). 
Temos, dessa forma, uma sequência de equações renormalizadas. 

Queremos agora obter uma relação análoga à Propriedade de Semigrupo. Note que por ()3.3|) é 
possível relacionar f n com a solução original: 

f n {x) = L a " + - +ai u (Z/ n+ '" +/3l x, L n ) . (3.6) 

Uma simples análise da equação (J3.6)) fornece a relação entre a dinâmica da sequência (f n ) quando 
n^ooeo comportamento assintótico de uq. Para isso, reescreva ()3.6|) como 

u o^^) = ^-J n (B nJ ^-), (3.7) 

onde A n e B n são chamados de prefatores e podem ser calculados por 

_ jjna n — (q„H hai) £j _ £n(3 n -([3 n -\ h/3i ) 

Assim, se os limites A n — > A, B n — > B, a n — * a, (3 n — > (3 e f n — > são obtidos quando n — ■> oo, 
podemos esperar que 

L na u Q {L n(5 x,L n ) = - An J n ( B n —^—-) — > A0(5x). 
O limite acima, por sua vez, estabelece a auto-similaridade do comportamento assintótico de uq. 
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3.2 O algoritmo 



A partir das definições preliminares, podemos sistematizar um método construtivo para a obtenção 
das sequências (a n ), (/3 n ), (A„), (B n ) e (/„), em que cada termo das sequências é calculado re- 
cursivamente. 

Sejam L, f e Gq tais que o problema de valor inicial (j3.1|) possui única solução u (x,t) para 
todo i 6 R e í > 1. O algoritmo iterativo consiste dos seguintes passos, para n = 1, 2, 3, . . ., 

1. Obtenha u n -\(x, L) evoluindo do dado inicial f n -\ do tempo t = 1 até í — L, de acordo 
com a equação ()3.4j) para u n -\. 

2. Calcule os expoentes a n e /9 n , conforme discutiremos mais adiante. 

3. Defina f n (x) = L^ u^L^x, L). 

4. Calcule os prefatores A n e B n por (|3.8|) . 

5. Calcule os parâmetros do operador G n a ser usado na próxima iteração, de acordo com a 
relação í)3.5j) . 

E o algoritmo numérico acima que usaremos para verificar as conjecturas já apresentadas ou até 
mesmo para se determinar os expoentes críticos sem que os mesmos sejam conhecidos a priori 
(como, por exemplo, na equação de Barenblatt). 

O expoente a pode ser sempre computado no passo Elpela fórmula a n = — ln^ w n _i(0, L) de forma 
que / n (0) = 1. A razão para isso é a suposta dinâmica auto-similar que esperamos encontrar: no 
regime assintótico, a solução u em x = se comporta como 



Como, por construção, u n (0, 1) = 1 para n — 1, 2, 3, . . ., e dada a relação entre u n no intervalo 
de tempo t G [1, L] e uq em t G [L n , L n+1 ] na equação ()3.3|) . temos que 



u (0,t) ~ c 



É 



— Q 



(3.9) 



no(0,£ n+1 ) 
«o(0,i>) 



f 
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e, por (J3.9J) . esperamos que, quando n — > oo, a n se aproxime de a. 

Enquanto o cálculo dos expoentes (a n ) pode ser sempre feito da forma descrita acima, a forma 
de se determinar os expoentes (f3 n ) depende do problema. Normalmente envolve uma relação de 
escala entre a n e (3 n de forma que uma certa parte do operador diferencial (escolhida a priori) 
permaneça invariante sob a renormalização da equação. 

Para ilustrar esta ideia, vamos rever brevemente o estudo da equação 

u t = u xx + Xu a u b x u c xx . 

A equação renormalizada é dada por 

ut = L l ~ 2 ^u xx + L^ b+2c ^ + ^ a ~ h ^ a ^Xu a ulu c xx . 

Se escolhermos (3 n = |, a parte linear da equação é mantida invariante e a dinâmica associada à 
equação 

u t = u xx (3.10) 

pode ser explorada. 

Por outro lado, se j3 n é escolhido de forma a satisfazer a relação 

l-{b + 2c)/3 n + {l-a-b-c)a n = 0, (3.11) 
então é o operador não-linear 

u t = Xu a u b x u c xx (3.12) 
que permanece invariante. Neste caso, é a dinâmica da equação ()3.12|) que pode ser explorada. 

Na linguagem do RG, a escolha de (3 n = \ concentra a atenção na dinâmica de ()3.1()j) e, dessa 
forma, é própria para a investigação de perturbações irrelevantes e marginais. A escolha 1)3.11)) 
leva o foco para ()3.12)) e permite a investigação dos casos de perturbações relevantes. 

Por último, a mudança de escala da variável espacial, no passo 01 pode ser feita de duas formas. 
A primeira foi utilizada por Chen e Goldenfeld [H] e consiste em modificar o tamanho da malha 
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Ax sem mudar os sítios discretos j — 0, 1, 2, . . ., de forma que depois da n-ésima iteração temos 
(Ax)' = L _ ^Ax e os pontos da nova malha estarão localizados em x = j(L _/3 Ax). A outra 
forma consiste em mudar os sítios discretos mantendo o tamanho da malha fixo, de forma que 
os sítios discretos, após a mudança de escalas, estarão em x = (L~^j)Ax. Neste caso, os valores 
da solução u nos pontos x = jAx da malha fixa têm que ser interpolados em cada iteração. 

A última técnica para se mudar a escala espacial referida acima foi proposta por Braga, Furtado 
e Isaia (veja a tese de doutorado de Isaia fH] e referências lá citadas). Utilizando este ponto de 
vista conseguimos um custo constante para cada iteração do algoritmo. Neste caso, temos um 
enorme ganho computacional, uma vez que, para se calcular u(-, L n ) a partir de u(-, 1) precisamos 
fazer n integrações no intervalo [1,L] a custo constante. 

3.3 Implementação 

Esta seção descreve algumas características da versão do algoritmo descrito acima que foi imple- 
mentada pelo autor desta dissertação. O algoritmo numérico foi implementado em linguagem C 
e a geração dos arquivos de entrada, bem como dos gráficos e tabelas a partir dos arquivos de 
saída, foi feita com o auxílio do software Maple©. Essa versão foi utilizada na obtenção dos 
resultados apresentados nos artigos |U 0] e no Capítulo |U desta dissertação. O leitor que esteja 
interessado apenas em conhecer o algoritmo em si, ou que queira implementar sua própria versão, 
pode dispensar a leitura desta seção. 

Equação e parâmetros 

O programa considera a equação 

u t = G(u,u x ,u xx ,x,t) = X (t p + 5f)[l + eH(-(u m ) xx )][l + ncos{ujx)]{u m ) xx + Xu a u b x u c xx , 

onde H é a função de Heaviside. Os parâmetros do programa são: L, x, P, õ, r, e, /x, cu, m, A, a,b e 
c, além da condição inicial, do número de iterações desejadas e da opção para mudança de escalas 
na variável espacial (veja este tópico mais adiante). Normalmente, escolhemos x = m = = 1 
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e e = O para o estudo dos Problemas ^ e El O parâmetro m é usado no estudo de equações 
que não consideramos nessa dissertação - portanto, em todas as simulações aqui mencionadas, 
assumimos a escolha de m = 1 - e o parâmetro e é usado para estudarmos a equação de Barenblatt 
na validação do método (veja a Seção l4~Tl no próximo capítulo). 



Discretização da equação 



Para a discretização da EDP, podemos utilizar qualquer esquema apropriado. Escolhemos um es- 
quema simples de diferenças finitas que combina o método de Euler para a discretização no tempo 
e, para a discretização no espaço, a fórmula padrão de três pontos para o operador Laplaceano 
e diferenças centradas para derivadas espaciais de primeira ordem. Dentro das limitações de 
estabilidade, o esquema resultante é de segunda ordem. 

Para discretizar o espaço-tempo, fazemos x = jAx e t = nAt, onde Ax é o tamanho da malha e 
Aí é o passo no tempo. Integramos a versão discreta da EDP para obtermos os dados (u(j Ax, L))j 
que aproximam a função u(x, L) nos sítios discretos x = jAx. Denotaremos por a aproximação 
para u(jAx, nAt). 

As derivadas serão aproximadas da seguinte forma: 

u n — 2u n + u n 

u xx (jAx,nAt) « J+1 \ 2 (3.13) 
u x (jAx,nAt) « J+ 2 1 (Ax) J " 1 . (3.14) 

Assim, temos o esquema de Euler explícito 

/ u n 

u] +1 = u] + AtG( u], 3+ *, A _>-\ 1 lK \., 1 ..iAv. nAt 



2(Ax) ' (Ax) 



Condição de estabilidade 

Reescrevemos a equação como 

u t = G(u, u x , u xx , x, t) = K(x, t, u xx )u xx + \u a u b x u c x 
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Uma condição suficiente de estabilidade será 1 : 

dG Ax dG > <9G <9G / Aí \ > 

dit xx 2 <9u x . ~~ <9u <9w xx \Ax 2 y ~~ 

Tome K', K" e R satisfazendo < A'' < A(x, í, u xx ) < A'" < oo. 
Supondo que A = 0, temos 

Aí 1 

como condição suficiente de estabilidade. 

Supondo que a solução u da equação satisfaça 2 \u\ < 1, |iia.| < 1 e (m^I < 1, temos as condições 
abaixo como condições suficientes para a estabilidade: 

A. < ^'-c|A|) 

Aí 1 

Ãã-2 - Ax 2 a|A| +2A // + 2c|A|' 

Condição de fronteira 

A vantagem do método explícito é que ele torna a velocidade de propagação da informação finita. 
Isso possibilita o estudo de problemas de valor inicial com dado inicial de suporte compacto 
adicionando-se mais um sítio discreto na variável espacial a cada novo passo de tempo. Dessa 
forma, a solução não é afetada por condições de contorno artificiais. 

Mudança de escala na variável espacial 

As duas formas de mudar a escala no passo El do algoritmo foram implementadas. A interpolação 

implementada foi a interpolação linear canónica. Uma desvantagem da interpolação é que ela não 

conserva a integral discreta da função interpolada. No estudo das equações lineares conservativas, 

em que o prefator depende apenas da integral do dado inicial, utilizamos a mudança de escala 
1 Veja o livro de Ames [T], Seção 2-18, p. 104. 

2 A limitação imposta sobre a solução u e suas derivadas é satisfeita por funções da forma <fi(x) = 
exp(— kx 2 /4:), k < 3, por exemplo. 
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sem interpolação para tornar o resultado mais fino. No Problema ^ com fi ^ 0, as oscilações 
do coeficiente de difusão têm seu período diminuído na mesma proporção que os intervalos da 
malha discreta no passo de mudança de escalas; neste caso, para enxergar o efeito da oscilação 
deste coeficiente com certa fidelidade, o recurso da interpolação também foi deixado de lado. Nas 
demais simulações utilizamos a interpolação. 

Ressaltamos que é possível fazer uma interpolação que conserve a integral discreta da função 
interpolada. Para isso deve-se adotar outro ponto de vista no momento da interpolação, que 
consiste em associar os valores às células e não aos pontos. Para o leitor interessado em se 
aprofundar nesta direção, indicamos o livro de LeVeque |21j . 

Fator de escala L 

A escolha do valor de L é puramente empírica. Em princípio qualquer fator de escala L > 1 
poderia ser usado. Para um valor de L muito grande, um elevado número de passos no tempo 
seria necessário para cada iteração do RG, e a malha espacial seria rapidamente contraída em 
poucas iterações devido à mudança de escala na variável espacial. Para um valor de L muito 
pequeno, seria necessário um número muito maior de iterações para se obter a acurácia desejada. 
Escolhemos um fator de escala L = 1,4 nas simulações em que utilizamos a interpolação no 
passo El escolhemos L entre 1,02 e 1,021 nos demais casos. 
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Capítulo 4 

Resultados Numéricos 



O objetivo deste capítulo é estudar os Problemas Q e El utilizando a implementação do RG 
numérico descrita no capítulo anterior. Primeiramente, procuraremos validar o método utilizando 
o mesmo no estudo de problemas de valor inicial em que soluções analíticas são conhecidas ou 
para os quais existem teoremas sobre a forma assintótica das soluções. Em seguida fazemos um 
estudo numérico da convergência do método, através de um refinamento de malha. Finalmente, 
usamos o algoritmo do RG para estudar os Problemas [T] e 121 do Capítulo [T] Cumpre notar que 
fazemos um estudo padronizado da equação do calor, em que verificamos: 

• a convergência do expoente a n para o valor teórico a = 1/2 com diferentes condições 
iniciais, de forma a confirmar a universalidade desse expoente crítico com respeito ao dado 
inicial; 

• a convergência do prefator A n para algum A > em todas a simulações feitas; 

• que o prefator A é de fato dado por (47r) -1 / 2 J r f(y)dy; 

• a convergência das funções de perfil nas normas L 1 e L°°; 

• que a função de perfil 0, obtida como lim^oo f n , é de fato a gaussiana 0*. 

Nos estudos subsequentes, alguns dos gráficos referentes aos itens acima serão omitidos por não 
trazerem grande novidade. As convergências de (B n ) e (J3 n ) são obtidas trivialmente uma vez 
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que o procedimento escolhe sempre o mesmo valor de (3 n para manter o termo linear invariante 
sob a renormalização da equação. 

A Seção E~Tl é dedicada à validação do método. Começamos pelo estudo padronizado da equação 
do calor, comprovando a universalidade do expoente a e da função de perfil <f), assim como a 
expressão para o prefator em termos da integral da condição inicial. Daí em diante, alguns 
gráficos análogos, que foram gerados em outras simulações numéricas, são omitidos por se tratar 
de repetições monótonas que acrescentam poucas informações, como a convergência da difer- 
ença relativa entre duas funções de perfil para zero ou a verificação de que a função de perfil 
é dada pela gaussiana. Estudamos então a equação de difusão não-linear e verificamos que as 
perturbações não-lineares classificadas como irrelevantes de fato não alteram o comportamento 
assintótico exceto pelo prefator A. A equação de Barenblatt é um exemplo bem ilustrativo: o 
expoente a não pode ser determinado de forma exata analiticamente nem formalmente, restando 
o recurso numérico como alternativa para se calcular este expoente. A validação neste caso é 
feita comparando-se o expoente calculado numericamente com aquele advindo de uma aproxi- 
mação de primeira ordem conhecida na literatura. Por último, escolhemos uma das condições 
iniciais e repetimos seu estudo com uma malha duas vezes mais fina. A semelhança visual entre 
os gráficos obtidos, exceto por uma ligeira suavização e maior precisão para a malha mais fina 
são os indicativos de convergência do método com relação ao refinamento da malha. 

Na Seção I4.2[ estudamos o problema de difusão com coeficiente periódico (Problema da In- 
trodução). Os resultados apresentados nessa seção são essencialmente aqueles feitos no artigo 
jl] e constituem-se principalmente da verificação da Conjectura 12.91 do Capítulo [21 O fato que 
traz mais novidade é a mudança da função de perfil devido à homogeneização do coeficiente de 
difusão. 

Na Seção IP1 é feito o estudo dos teoremas e das conjecturas da Seção l2~!T2l do Capítulo que 
tratam do problema de difusão com coeficiente dependente do tempo (Problema|2]da Introdução). 
Esse estudo pode ser encontrado no artigo em elaboração [3]. As parte mais interessantes, na 
opinião do autor, são a ampliação da classe de perturbações não-lineares irrelevantes parap > 0, 
que passa a incluir todas aquelas que eram marginais nas outras equações, como u 3 , uu x e u xx , 
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e a transição de fase entre os regimes assintóticos linear e nao-linear para uma forma particular 
da equação com expoente não-inteiro e solução positiva. 

Como dito acima, os gráficos que trazem informação repetitiva são omitidos nas Secões l4~2*l e l4.3l 

Em todo este capítulo, o símbolo <ft vai denotar a função f n obtida no último passo da iteração 
do algoritmo, após terem sido obtidas as convergências desejadas. 

Chamamos a atenção para uma diferença na definição das funções de perfil. A versão numérica 
que implementamos (descrita no Capítulo EJ), calcula a n de maneira a forçar que / n (0) = 1 para 
n > 1. Por este motivo, será conveniente redefinir 

Isso implica na necessidade das constantes ^Jl/Ait, \J {p + l)/4ir e ^l/Aira, respectivamente, 
multiplicando o prefator A dos teoremas e conjecturas apresentados no Capítulo Ele mencionados 
na Introdução. Em particular, teremos as seguintes relações, que serão usadas mais adiante: 

-ln& = x 2 /4 (4.1) 
-ln0 p = (p + l)-(x 2 /4) (4.2) 
ln = --(x 2 /4). (4.3) 



4.1 Validação do método 
4.1.1 Equação do calor 

Vimos no Capítulo |2 que o comportamento assintótico das soluções da equação do calor com 
condição inicial u(-, 0) = / £ Ba é dado por 

u (a;,í)«^(£), (4.4) 

com a = (3 = 1/2, <p = (fi* e A = (47r) -1 / 2 J R f(y)dy. Segundo os resultados teóricos já apresen- 
tados, a, (3 e são universais com respeito à condição inicial, desde que a mesma esteja, por 
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exemplo, em B4. Por outro lado, o prefator A depende apenas do que chamamos de massa do 
dado inicial /, dada por Lf(y)dy. 

As três condições iniciais utilizadas, fi, f% e fa, são exibidas nas Figuras l4~Tj 14.21 e 14. 3| respecti- 
vamente. 



/ 0.8- 




/ 0.6- 




/ 0.4- 




/ 0.2- 

' ~ 1 1 1 1 1 l 1 1 1 


1 1 1 1 | 1 1 1 1 | ! ! 



Figura 4.1: Gráfico da condição inicial fi. 

A Figura l4~4l mostra a convergência do expoente a n , calculado na n-ésima iteração do algoritmo, 
para o valor teórico a = 1/2, com diferentes condições iniciais. Isso ilustra a universalidade de a 
com respeito a condições iniciais suficientemente localizadas (por exemplo, de suporte compacto). 

A Figura l4~51 mostra a convergência dos prefatores A n para os valores teóricos. 

A Figura 1431 mostra uma curva parametrizada cujas componentes são (x 2 /4, — ln<j)(x)). O fato 
de todo ponto desta curva estar em uma reta de inclinação igual a 1 significa que = 0*. 

Finalmente, na Figura l4~7l plotamos a diferença relativa entre duas funções de perfil consecutivas, 
Wfn — /ra-i||/||/ra||; como uma função do número n de iterações do RG. As normas que utilizamos 
em todas as simulações são L 1 e L°°. O gráfico exibido na Figura I4~7I foi gerado com a norma 
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Figura 4.2: Gráfico da condição inicial fa- 

L 1 (M). A convergência dessa diferença para zero é forte indicativo da convergência de (f n ) 
e, consequentemente, de um comportamento assintótico auto-similar. O cálculo da diferença 
relativa também pode ser utilizado como um critério de parada das iterações. 
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O 100 200 300 400 

Figura 4.5: Curvas A n x n para a equação do calor com diferentes condições iniciais. As linhas 
pontilhadas correspondem aos valores teóricos. Estes valores, dados por (4ir)~ lj/2 f R f(y)dy, são 
1,000, 1,133 e 0,8901. 




Figura 4.6: Curvas — ln0 x x 2 /4 para a equação do calor com diferentes condições iniciais (as 
curvas se sobrepõem). 
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Figura 4.7: Curvas (||/ n — f n -i\ \/\ \f n \ |) x n para a equação do calor com diferentes condições 



iniciais. 
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4.1.2 Equação de difusão não- linear 

O Teorema 12.71 diz que soluções da equação 

u t = u xx + Xu a u b x u c xx , 
com condição inicial u(-, 0) = / G £>4 e a + 26 + 3c — 3 > 0, se comportam como 

u(z,í)«^(£), (4.5) 

com a — (3 — 1/2, (p = </>*, e A é uma constante positiva que pode depender de A, a, b, c e c. 

Fizemos simulações numéricas em que variamos a condição inicial, o parâmetro A e os expoentes 
a, b e c de maneira a satisfazer à condição a + 2b + 3c — 3 > 0. As simulações são descritas pela 
Tabela WÃ\ e a universalidade de a é ilustrada na Figura l4~8l 

Para todas as simulações, verificamos que <fi = 0*, que o prefator converge para algum número 
positivo e que a norma relativa da diferença entre duas funções de perfil consecutivas converge 
para zero. 



n 


A 


a 


b 


c 


1 


.10 


3 


1 





2 


.20 


3 


1 





3 


-.10 


1 


1 


1 


4 


.10 





1 


1 


5 


.15 


3 


1 





6 


-.20 


1 


1 


1 


7 


.30 


1 





1 



Tabela 4.1: Descrição das simulações feitas no estudo das classes de universalidade do compor- 
tamento assintótico das soluções da equação de difusão não-linear, numeradas de 1 a 7. 
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Figura 4.8: Gráfico de a n x n para as simulações descritas na Tabela l4~T| ilustrando a universal- 
idade do expoente crítico a para perturbações irrelevantes da equação de difusão não-linear. 
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4.1.3 Equação de Barenblatt 



Estudamos também a equação de Barenblatt 

ut=[l + eH(-u t )}u xx , (4.6) 

onde H(v) é a função de Heaviside com salto em v — 0. Resultados teóricos 1 fornecem uma 
aproximação de primeira ordem para o expoente a dada por 

Simulações foram feitas para diferentes valores de e. Os valores de a calculados numericamente 
são comparados com o resultado teórico acima na Figura PPfl 

Em todas as simulações, verificamos que o prefator A n converge a um número positivo e a norma 
relativa da diferença de duas funções de perfil consecutivas converge a zero. 

O estudo numérico desta equação também ilustra uma característica interessante do RG numérico: 
é possível calcular os expoentes críticos ainda que os mesmos não sejam conhecidos de forma exa- 
ta. 



1 Ver ^H] e referências lá citadas. 
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0.55 



0.5 



0.45 



0.4 -I ;"' 

-0.4 -0.2 0.2 0.4 

Figura 4.9: Curvas a x e para diferentes simulações com a equação de Barenblatt, onde a é 
calculado numericamente e a reta pontilhada é dada pela aproximação de primeira ordem. 
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4.1.4 Refinamento da malha 



Para a equação do calor, fizemos um teste de refinamento de malha utilizando a condição inicial 
/2 e uma outra idêntica à /2 porém com uma malha duas vezes mais fina. Os gráficos de a n e A n 
são exibidos nas Figuras 14.101 e 14.111 respectivamente. Na simulação com a malha refinada foi 
verificado que = 0* e que a norma relativa da diferença de duas funções de perfil consecutivas 
converge a zero. 

0.53 




400 



Figura 4.10: Curva a n x n para a equação do calor com condição inicial uma curva foi gerada 
com a malha duas vezes mais fina que a outra. 
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Figura 4.11: Curva A n xn para a equação do calor com condição inicial f'2] uma curva foi gerada 
com a malha duas vezes mais fina que a outra. 
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4.2 Difusão com coeficiente periódico 



As simulações discutidas nesta seção visam a corroborar a Conjectura 12.91 

Nas simulações realizadas, observamos a convergência de a n para a = 1/2. Para verificar a 
universalidade deste fato, utilizamos diferentes condições iniciais (Figura 14.12)) . perturbações 
não-lineares ÍFigura l4.14|) e valores de \i ÍFigura l4.13J) . 



Em todas as simulações, a convergência da norma relativa da diferença entre funções de perfil 
consecutivas convergiu para zero e o prefator A n convergiu para um número positivo. 



1 

0.8 
0.6 
0.4 
0.2 
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Figura 4.12: Curvas a n x n para o problema de difusão com coeficiente periódico e diferentes 
condições iniciais. 
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Figura 4.13: Curvas a n x n para o problema de difusão com coeficiente periódico e diferentes 
valores de \i. 

0.9- 
0.8: 
0.7: 
0.6 :v 

0.5:^^SS^ 

0.4: / 

0.3:/ 

0.2: 

°- 1 50 100 150 200 250 300 350 400 

Figura 4.14: Curvas a n x n para o problema de difusão com coeficiente periódico e diferentes 
perturbações não-lineares irrelevantes. 
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4.2.1 Homogeneização do coeficiente de difusão 

A principal diferença desta equação em relação à equação do calor, ou à equação de difusão 
não-linear, é a mudança da função de perfil causada pelo termo [ig(x). Por isso, plotamos as 
curvas — \n(<f>(x)) x x 2 /4 para diversos experimentos. O fato de estas curvas serem retas com 
inclinações a~ l significa que <fi = cf) a para cada uma delas - veja (JUBJ). A Figura E. 151 mostra essa 
curva para fi = 0, 8 e g — cos, juntamente com o valor de o = 0, 6 esperado teoricamente. 




Figura 4.15: Curva parametrizada (— ln0, x 2 /4), cujo traço é uma reta com inclinação a 
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4.3 Difusão dependente do tempo 



Nesta seção vamos estudar os resultados apresentados na Seção 12.3.21 do Capítulo [21 Esses 
resultados consideram o problema de difusão com coeficiente dependente do tempo. O caso linear 
é coberto pelo Teorema 12.111 enquanto para o caso não-linear foi formulada a Conjectura 12.121 
Começamos pelo caso linear. 

4.3.1 Caso linear 

Considere a equação 



onde r < p, com condição inicial u(-,l) — f E B4. O Teorema 12.111 diz que as soluções deste 
problema se comportam assintoticamente como 



Para estudar numericamente o resultado acima, realizamos um conjunto de 18 simulações, que 
correspondem à combinação de três condições iniciais diferentes com coeficientes de difusão de 
seis formas diferentes, dadas por: í 1 / 2 + í 1//4 ; í 1//2 + 1; í 1 / 2 ; t + í 3//4 ; t + 1; e t. 

A Figura E.16I mostra a convergência dos expoentes a para o valor (p + l)/2 e sua universali- 
dade com respeito ao termo ôt r e à condição inicial /. Os valores de p utilizados são 1/2 e 1, 
determinados pelos coeficientes K(t) listados no parágrafo anterior. 

O prefator A é também universal com respeito a detalhes nas condições iniciais e depende apenas 
da integral das mesmas. A Figura i!. 171 mostra a convergência do prefator para os valores teóricos 




coma = f3=(p + l)/2, = P e A = [(p + 1)/A^ 2 J R f(y)dy. 



A=y/(p + l)/47rJ R f(y)dy. 



Por ultimo, a Figura EHH1 revela que a função de perfil é, de fato, dada por = <f> p . 
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10 20 30 40 50 60 



Figura 4.16: Curvas a n x n para a equação linear com coeficiente dependente do tempo para 
diferentes coeficientes e condições inicias. Os valores teóricos que correspondem a p = 1/2 e 
p—l são, respectivamente, a = 0, 75 e a = 1, 0. 




Figura 4.17: Curvas A n xn para a equação com coeficiente dependente do tempo para diferentes 
coeficientes e condições inicias. As linhas pontilhadas correspondem aos valores teóricos. Estes 
valores, dados por [(p + 1)J^k\ 1/2 J R f(y)dy, são 1,090, 1,225, 1,259, 1,388, 1,414 e 1,602. 
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Figura 4.18: Curvas parametrizadas — ln0 x x 2 /4 para a equação com coeficiente dependente 
do tempo com diferentes coeficientes e condições inicias. A inclinação da reta deve ser dada por 
p + 1. Cada curva visível no gráfico é na verdade a superposição de nove curvas geradas pelas 
simulações descritas acima. As linhas pontilhadas e o ponto marcado (l,p + 1) correspondem 
aos valores teóricos. 
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4.3.2 Caso não- linear 

Consideramos a equação 

u t = (t p + 5t r )u xx + Xu a u b x u c xx . 

Todas as perturbações não-lineares que foram classificadas como irrelevantes para a equação 
de difusão não-linear, ou seja, aquelas que satisfazem cLf = a + 2b + 3c — 3 > 0, ainda são 
irrelevantes para a equação com coeficiente de difusão dependente do tempo para p > 0, pois 
satisfazem rjp = a + 2b + 3c — (3 + p)/(l + p) > 0. 

No caso em que p > 0, as perturbações que antes eram classificadas como marginais, agora 
passam a ser também irrelevantes. Essas perturbações podem ser F = u 3 , F = uu x ou F = u xx . 

A Conjectura 12.121 diz que o comportamento assintótico das soluções da equação acima ainda é 
o mesmo dado pela equação linear tratada pelo Teorema 12.111 exceto pelo fato de que o prefator 
A é uma constante positiva que depende da condição inicial, do coeficiente de difusão e da 
perturbação não-linear, desde que a perturbação não-linear F satisfaça r/p > 0. 

Os gráficos relativos às simulações feitas nesta seção são omitidos por serem todos análogos aos 
que já foram apresentados em seções anteriores. Por esta razão vamos nos limitar a exibir a 
Tabela E~2"l que descreve as simulações realizadas visando a confirmar a afirmação do parágrafo 
anterior. 
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n 


P 


A 


a 


b 


c 


1 


.50 


.10 


1 


1 





2 


.50 


.10 


3 








3 


.50 


.10 








1 


4 


.25 


.10 


1 


1 





5 


.25 


.10 


3 








6 


.25 


.10 








1 


7 


.50 


.10 


3 


1 





8 


.50 


.20 


3 


1 





9 


.50 


-.10 


1 


1 


1 


10 


.50 


.10 





1 


1 


11 


.50 


.15 


3 


1 





12 


.50 


-.20 


1 


1 


1 


13 


.50 


.30 


1 





1 



Tabela 4.2: Descrição das simulações feitas no estudo das classes de universalidade do compor- 
tamento assintótico das soluções da equação de difusão não-linear com coeficiente de difusão 
dependente do tempo, numeradas de 1 a 13. 
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4.3.3 Transição de fase 



Vamos estudar numericamente a equação 



Uf tPu xx 



u 



a 



(4.7) 



onde a e p sao parâmetros reais tais que a > 1 e p > 0. 

Dado p > 0, o expoente a que corresponde a uma perturbação marginal é 2 



No plano (p, a), a curva ()4.8)) é formalmente a interface entre a região relevante a > (í>+3)(p+l) _1 
e a região irrelevante a < {p + 3)(p + l) -1 . O estudo da função de perfil, do expoente crítico a e 
dos prefatores, indica que ()4.8)1 é de fato a curva crítica para a transição de fase entre os estados 
relevante e irrelevante. 

A Figura l4*. 19l mostra o espaço de parâmetros (p, a) com a curva crítica (ou separatriz) juntamente 
com as simulações que fizemos neste trabalho. A curva a c (p) = (p + 3)(p + 1) _1 corresponde 
aos valores de (p, a) que tornam a equação ()4.7|) invariante à renormalização com expoentes 



Para os valores de (p, a) acima da curva crítica, o operador não-linear tende a zero sob a ação de 
sucessivas renormalizações com (3 = a = (p+ 1)/2, enquanto a parte linear permanece invariante. 
Neste caso, o regime assintótico das soluções é essencialmente da equação linear. 

Para os valores de (p, a) abaixo da curva crítica, a equação (|4.7jl permanece invariante para 
(3 = (p+ l)/2 and a — l/(a — í). Neste caso o expoente de decaimento é determinado pela parte 
não-linear do operador, enquanto a difusão é devida ao termo linear, e a função de perfil depende 
de ambos. 

A Figura E.20I mostra o gráfico da função 



(3 = (p+l)/2 e a 



(p + l)/2 = l/(o-l). 




(4.9) 



2 Faça i] = cm l|1.10H ■ 
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Figura 4.19: Espaço de parâmetros (p, a) para a equação (|4.7jl . A curva crítica separa os regimes 
linear e não-linear da equação. Os quadrados representam os parâmetros para os quais foram 
feias as simulações numéricas. 

Olhando para as curvas de nível da superfície gerada, é possível perceber que as singularidades 
da superfície (p, a, a) formam uma curva que corresponde aos parâmetros peada curva crítica. 
E possível também perceber que de um lado dessa curva o expoente a é determinado por p 
enquanto do outro lado é determinado por a. 

Com o objetivo de verificar a relação (j4.9|) . na Figura 11.211 plotamos a interseção da superfície 
()4.9|) com o plano p = 1/2, juntamente com os expoentes a calculados pelas nossas simulações. 
Além disso, a Tabela l4~31 compara os valores teóricos com os valores computados de a como uma 
função de a para p = 1/2 fixo. A Figura E. 221 e a Tabela l4~4l são análogos para a = 7/3 fixo e a 
como uma função p. 

As Figuras 14.231 e l4~.24l ilustram a diferença na convergência de (a n ) e (A n ). Fazendo p — 1/2 
e escolhendo três valores diferentes para a (um crítico, um subcrítico e supercrítico), podemos 
investigar a diferença dos regimes, no que diz respeito os principais aspectos qualitativos de seu 
comportamento assintótico. 

Na Figura H. 231 podemos ver a diferença no cálculo do expoente a. No caso supercrítico (termo 
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Superfície alfa(p,a) 




Figura 4.20: Superfície a(p,a) dada por (|4.9|) . As singularidades da superfície correspondem à 
curva crítica sobre a qual ocorre a transição de fase. 

não-linear irrelevante), o expoente de decaimento é determinado por a = [(p + l)/2] = 3/4, 
enquanto que no caso subcrítico este expoente assume um valor maior. Além disso, quando 
usamos parâmetros críticos, a convergência (a n ) para a = [{p + l)/2] = 3/4 é muito mais lenta 
que nos outros casos (apesar de não aparecer na figura, após um número muito maior de iterações 
podemos ver que a convergência ocorre de fato). 

Na Figura II. 241 vemos a convergência do prefator nos três regimes. Para parâmetros na região 
relevante e irrelevante, o prefator converge rapidamente a um número positivo, enquanto no caso 
crítico o prefator converge para zero quando n — > oo. Este fato é devido ao aparecimento de uma 
correção logarítmica ao comportamento auto-similar, que será discutido logo adiante. 

Na Figura E.25I ilustramos a universalidade do prefator com respeito à condição inicial no caso 
de perturbações relevantes (a subcrítico). Fixamos p — 1/2 e escolhemos três valores distintos 
de a satisfazendo a < a c (p). O prefator obtido como lim n ^oo A n depende apenas da escolha de a 
e não das condições iniciais. 

A Figura E.26I mostra a não-universalidade do prefator no caso da perturbação irrelevante. Para 
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Figura 4.21: Gráfico de a x a para p — 1/2 fixo. Os quadrados representam os valores calculados 
numericamente e a linha contínua é dada pela equação (|4.9|) . 

p = 1/2 fixo, escolhemos três valores de a > a c (p) e três condições iniciais distintas. O prefator 
obtido é influenciado tanto pelas escolha de a quanto da condição inicial. 
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a 


a teórico 


a calculado 


2.133333 


.882353 


.882353 


2.183333 


.845070 


.845071 


2.233333 


.810811 


.810811 


2.283333 


.779221 


.779221 


2.333333 


.750000 


.750553 


2.383333 


.750000 


.749991 


2.433333 


.750000 


.749991 


2.483333 


.750000 


.749991 


2.533333 


.750000 


.749991 



Tabela 4.3: Valores de a como uma função de a para p — 1/2 fixo. O valor calculado numerica- 
mente é comparado com o valor dado pela equação ()4.9|) 



p 


a teórico 


a calculado 


.300000 


.750000 


.750000 


.350000 


.750000 


.750000 


.400000 


.750000 


.750000 


.450000 


.750000 


.750001 


.500000 


.750000 


.750553 


.550000 


.775000 


.774994 


.600000 


.800000 


.799996 


.650000 


.825000 


.824987 


.700000 


.850000 


.849962 



Tabela 4.4: Valores de a como uma função de p para a = 7/3 fixo. O valor calculado numerica- 
mente é comparado com o valor dado pela equação ()4.9|) 
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074 0.3 0.4 0.5 0.6 0.7 

P 

Figura 4.22: Gráfico de a xp para a = 7/3 fixo. Os quadrados representam os valores calculados 
numericamente e a linha contínua é dada pela equação ()4.9j) . 




Figura 4.23: Curvas a n x n para três regimes distintos. Fixado o parâmetro p = 1/2, cada curva 
foi gerada com um valor de a < 7/3, a = 7/3 e a > 7/3, correspondendo aos regimes não-linear, 
crítico e linear, respectivamente. 
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Curvas A n x n 




200 400 n 600 800 ^1000 



Figura 4.24: Curvas A n xn para três regimes distintos. Fixado o parâmetro p — 1/2, cada curva 
foi gerada com um valor de a < 7/3, a = 7/3 e a > 7/3, correspondendo aos regimes não-linear, 
crítico e linear, respectivamente. 
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Curvas A x n para o regime nao-linear 
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Figura 4.25: Curvas A n xn para o regime não-linear, com três valores distintos de a satisfazendo 
a < 7/3 e três condições iniciais distintas, fixado p — 1/2. O prefator dado por A = lim n ^oo A n 
depende apenas de a. 
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Curvas A x n para o regime linear 
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Figura 4.26: Curvas A n x n para o regime linear, com três valores distintos de a satisfazendo 
a > 7/3 e três condições iniciais distintas, fixado p = 1/2. O prefator dado por A = lim n ^oo A n 
depende do parâmetro a e da condição inicial utilizada. 



79 



4.3.4 Correções logarítmicas sobre a curva crítica 



Para a equação (J4.7|) com parâmetros sobre a curva crítica, o comportamento assintótico auto- 
similar esperado sofre uma correção logarítmica: 

u(x, t) ~ -7— — r— <p { ~n x J quando t — > oo, (4-10) 



(ílní) a yn 3 

onde (3 — (p + l)/2 e a = (p + l)/2 = l/(a — 1). Além disso, resultados teóricos dizem que 
o prefator A não depende da condição inicial, mas apenas da escolha dos parâmetros (p, a). 
Lembrando que o algoritmo faz com que /„(0) = 1 no processo construtivo, combinando (j4.1(J|) 
com (|3.7|) obtemos 

ln A n m ln — — alnn. (4-11) 

Essa relação diz que o gráfico de lnA n x lnn se aproxima, para n suficientemente grande, de uma 
linha reta tendo inclinação dada por —a que intercepta o eixo vertical no ponto \nA/L a . Neste 
caso, a universalidade de A traduz-se no fato de que quaisquer condições iniciais dentro de uma 
classe apropriada dão origem a curvas que são assintóticas à mesma reta e não a retas paralelas. 

As afirmações do parágrafo acima foram verificadas para p = 0, p = 1/2, e p = 1. A Figura E. 2 71 
mostra o gráfico de lnv4 n x Inn para p — 1/2 e três condições iniciais diferentes. A inclinação da 
reta assintótica é —a, sendo utilizado o valor de a calculado pelo RG. 
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Curvas log A x log n 



1 




Figura 4.27: Curvas lnA n x lnn parap = 1/2 e a = 7/3, com três condições iniciais distintas. A 
reta pontilhada passa pelo ponto (ln.A n ,lnn) da última iteração do algoritmo e tem inclinação 
—a, considerando o expoente calculado na última iteração. 
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Capítulo 5 

Teoria de Homogeneização 



Neste capítulo vamos apresentar a teoria de homogeneização à qual nos referimos nos capítulos 
anteriores. A principal ideia na teoria de homogeneização é a de tentar descrever um meio não- 
homogêneo a partir de seu comportamento efetivo. Essa tentativa baseia-se na hipótese de que 
existe um meio homogéneo cujas propriedades são próximas daquelas do meio real quando me- 
didas em grandes escalas. O processo de "aproximar pela média", chamado homogeneização, 
substitui uma estrutura complicada em escala microscópica por uma estrutura homogénea assin- 
toticamente equivalente. A teoria que vamos apresentar neste capítulo diz de uma forma matem- 
aticamente precisa que a substituição da equação original por outra com coeficiente constante é 
uma aproximação válida em relação a um certo limite. 

Escrevemos esse capítulo com dois objetivos. O primeiro é chegar ao Teorema 15.1 li e ao Corolá- 
rio que justificam a substituição do coeficiente não-homogêneo por um homogéneo no estudo 
do Problema ^ da Introdução. O segundo é apresentar ao leitor um texto auto-contido sobre a 
teoria de homogeneização, para equações elípticas e parabólicas lineares, em espaços de Sobolev. 
O conteúdo deste capítulo se baseia no artigo de Papanicolaou [21] e no livro de Jikov et al. |2Uj . 

Vamos fazer uso indiscriminado dos resultado básicos da teoria de medida e integração em subcon- 
juntos de M. m . Serão utilizados também alguns elementos de análise funcional como os conceitos 
de espaços de Hilbert e de Banach, ou o Lema de Riez. Contudo, não é necessário que o leitor 
tenha feito um curso de análise funcional desde que tenha algum livro introdutório para consultar. 
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Na Seção 15.11 vamos apresentar as ideias básicas da técnica de homogeneização e estudar um 
caso simples em que não são necessários muitos resultados matemáticos mas que, por sua vez, 
contribuem para uma visão simples e clara do problema de homogeneização. 

Na Seção 15.21 vamos definir alguns espaços de Sobolev e suas principais propriedades, trabalhar 
com seus duais e estudar a equação elíptica linear com condição de contorno periódica. Isso feito, 
introduzimos um problema auxiliar, que prepara o terreno para o Teorema 15.111 em que uma 
equação parabólica mais geral do que o Problema ^ linear é estudada. 

5.1 Estudo de um problema unidimensional 

Nesta seção pretendemos apresentar o estudo de um caso bastante simples com o intuito de dar 
um enfoque físico e intuitivo da técnica de homogeneização. Por esta razão, optamos por deixar 
em segundo plano as dificuldades matemáticas (utilizando para isso o recurso de acrescentar 
mais hipóteses do que o necessário aos problemas), e demos preferência à tentativa de justificar 
e ilustrar a idéia de homogeneizar um meio heterogéneo. 

O problema elíptico linear unidimensional 

Considere a equação abaixo, que modela o estado de equilíbrio de um condutor (de calor) unidi- 
mensional ao qual aplicamos um gradiente de temperatura V: 



onde N > é a metade do tamanho do condutor que estamos considerando e D(x) é uma 
função periódica (de período 1), limitada e satisfaz mí x ^D(x) > 0. A função D(x) representa 
a condutividade térmica no ponto x e V representa o gradiente de temperatura aplicado ao 
condutor. 

O fluxo em um ponto x é dado por _D(x)^— ; o fluxo médio, por unidade de gradiente de 




±-L = 0, xe [-N, N] 
ix 

VN, u(+N) = +VN 
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temperatura, será denotado por Dn e dado pela fórmula 



V \2N J_ N y ' dx 

No caso unidimensional, que estamos tratando neste momento, o fluxo D(x)u x (x) é constante, 
logo a integral acima é desnecessária. Entretanto, essa representação tem o significado físico de 
estarmos tomando a média do fluxo. Em dimensões maiores, ao invés de um fluxo constante 
temos um fluxo com divergente nulo. 

Queremos determinar a condutividade efetiva D*, que é definida pela existência do limite 

D* = lim D N . 

AT— >oo 

Ressaltando novamente que estamos tratando do caso unidimensional, temos uma fórmula ex- 
plícita para D N . Cada equação abaixo é consequência imediata da equação que a precede: 

D N = -D(x)u x (x) Vxe[-N,N] 
u x = VD N /D(x) 

+N n+N J 

2VN= I u x dx = VD N / -rr^dx 

Dlx) 



N J -N 



1 r+ N i 



Como D é periódica, de período 1, por exemplo, a última igualdade implica em 

1 , \ /— i\-i 



Uma interpretação física para a fórmula acima é que a resistividade efetiva (D*) -1 é a resistividade 
média (D" 1 ) em um sistema de pequenos resistores em série. Observamos que o mesmo raciocínio 
pode se aplicar ao caso em que a conditividade D não é periódica ou até mesmo quando é um 
campo estocástico, exceto pela integral desta última fórmula. 



Uma forma equivalente de ver o mesmo problema 



O mesmo problema pode ser visto sem que se aumente o tamanho do condutor no limite N — » oo, 
como fizemos acima. Podemos ajustar nossa escala de medida ao tamanho do problema e pensar 
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que, nessa escala, a heterogeneidade é que passa a ter uma oscilação cada vez mais rápida. Isso 
é feito através de uma mudança de variáveis. Primeiramente fazemos 

1 

e definimos 

D e (x) = D(x/e) e u e (x) = eu{x/e). 
A nova equação a ser considerada será 

" [^(4 = *e [-1,1] 



e o coeficiente efetivo dado por 

1 r 1 [àu. e (v\~\ 

dx. 

áx 

Com essa abordagem obtemos o mesmo valor para D* 



±1 






du e (x) 


f D<{x) 




dx 



O problema que vamos considerar em seguida ilustra bem a ideia de que a aproximação do meio 
heterogéneo pelo meio homogéneo é uma "boa" aproximação. 

Um problema de homogeneização 

Considere D e C 1 (1R) positiva e periódica de período 1. Para cada e > defina D € (x) = D(x/e) 
e seja u e a solução de 

/ - [D*{x)u%{x)\ x = /, X6[-1,1] 
\n e (±I) = ' 

onde/ eC^I-M]). 

Seja u° a solução do mesmo problema substituindo-se D e por 

" = (I m*T = {Drl >~ 1 ■ (5 - 2) 

O restante desta seção é dedicado ao propósito de oferecer uma prova simples e elegante da 
seguinte proposição. 
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Proposição 5.1 Com u e e u definidas acima, vale o limite 



u e (x) — >u°(x) (5.3) 



uniforme em [—1,1]. 



A equação unidimensional (jõ.ljl pode ser resolvida explicitamente. Integrando a equação de a 
x obtemos 

D%x)u%{x) = [£> £ (0K(0)] - [ X f(t)dt. 

Jo 

Definindo F(x) = f(t)dt, temos 

Integrando novamente, obtém-se 

u\x) = u e (0) + [D e (0)<(0)] / ttVt^ - / 7TTT dí - 



, £>«(*) y D*(t) 
Definimos Ci(e) = [D e (0)tt^(0)] e c 2 (e) = w e (0), donde 

u\x) =c 2 (e) + Cl (e) f X -L-dt- P Py^dt. (5.4) 



Jo Z>(í) y D«(í) 

A condição de contorno é u e (±l) = 0. Subtraindo ou somando as equações u e (— 1) = e 
u e (+l) = obtemos, respectivamente, 



! y y_ x ^(t) 



. . 1 / f +1 sgn(í) N , , , sgn(í) , 



onde sgn(0) = e sgn(x) = x/\x\ para i^O. 

Como vP satisfaz a mesma equação com D e substituído por D*, fazemos a mesma substituição 
em ()5.4|) para concluir que u° será dado por 

u°(x) = c 2 (0) + Cl (0)-^ - jT F(f)dí, (5.5) 
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com 



1 



c 2 (0) = 2D^J sga(t)F(t)dt 



Valem as seguintes convergências uniformes em [—1,1]: 



1 -dt— e f X iMLdt—^^- I F(t)dt. (5.6) 



/ D'(t) e^o D* J Q D £ (t) e-o D* J 

Para provar essas convergências, vamos utilizar uma propriedade geral de funções periódicas, 
dada pelo lema abaixo. Resultados análogos valem sob condições bem mais gerais, mas com a 
hipótese de suavidade a demonstração se torna particularmente simples. 

Lema 5.2 (Propriedade do Valor Médio) Seja F : R x [a, b] — > IR de classe C 1 , tal que 
F(y + 1,2;) = F(y,x) para todos x G [a,b], y G R. Então, dados a',b' G [a, 6], a' < i>a/e o 
seguinte limite 

lim / F(x/e,x)dx = / F(x)dx 

40 y a , y a , 



uniforme em relação a a! e b' , onde 



F(x) = / F(y,x)dy. 
Jo 

Prova: Considere G:lx[a,&]^R dada por 

rv 

G(y,x) = / [F(z,x)-F(x)]dz. 
Jo 

Temos G de classe C 1 , periódica de período 1 em y, pois é a integral de uma função periódica da 
qual subtraímos a média; portanto G é limitada. Temos que 

d 1 G(y,x) = F(y,x)-F(x). 

Segue que 

-6' /•&' 

[F(x/e, x) — F(x)]dx = / d\G{x/e,x)dx 

Ja' 
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V/e 



diG(y, ey)edy 



a'/e 



b'/e çV/e 

[d 1 G{y,ey) + ed 2 G{y,ey)]dy- / d 2 G(y,ey)edy 

a'/e Ja'/e 
V/e d rV/e 

—G{y,ey)dy- / d 2 G{y,ey)edy 

a'/e a V Ja'/e 



G(b'/e, b') - G(a'/e, a') 
Aplicando a desigualdade triangular, temos 

F(x/e, x)dx 
donde temos o limite uniforme desejado. 



d 2 G(x/e, x)dx 







/ F{x)dx 


< e 


la' 





2 max \G(y, x)\ + (b — a) max \d 2 G(y,x) 
[o,i]x[a,6]' vy ' Jl v ; [o,i]xm' vy ' ' 



□ 



Pela definição de D* dada por (|5.2jl . a Propriedade do Valor Médio tem como consequência 
imediata que vale o limite 

[ ik) gm ^ [ g(t)it 

uniforme em [—1,1], onde g é qualquer função em C 1 ([— 1, 1]). 



Em particular, temos que valem os limites (|5.6[) . uniformemente em [—1,1]. Também valem os 
nites Ci(e) — 

€ — I 

-1,0) e (0,1]. 



limites ci(e) — > ci(0) e c 2 (e) — > c 2 (0), separando-se as integrais de c 2 (e) e c 2 (0) nos subintervalos 

e— >0 e— »0 



Comparando (|5.4|) e (|5.5|) . temos o limite ()5.3|) uniforme em [—1,1] como consequência dos limites 
mencionados acima. 



5.2 Equações diferenciais com coeficientes periódicos 
5.2.1 Espaços de Sobolev 

Nesta seção Q denotará um conjunto aberto limitado em M m . Os conjuntos de funções reais 
L a (Q), 1 < a < +oo são definidos da maneira usual. O conjunto C^°(Q) é formado pelas 
funções reais infinitamente diferenciáveis com suporte compacto em Q. 
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Definição 5.3 (Gradiente Fraco) Uma função v : Q — > M. m , dada por v = (v\, . . . ,v m ), V{ G 
L l (Q) é chamada gradiente da função u G L l (Q) se 

d f 

u(x)— — <f>{x)dx = — / vAx)(j){x)dx, V0 G C^(Q), i — l,...,m. (5.7) 
OXi J Q 

Denotamos o gradiente por Vw ou J^u. 

Observe que o gradiente, quando existe, é único q.t.p. e que V é linear, isto é, V(«i + Au 2 ) = 
Vui + ÀVw2, sempre que V«i e V«2 existirem. Por exemplo, se u, v : (—1, 1) — > R são dadas por 
u(x) = \x\ q.t.p. e v(x) = x/\x\ q.t.p. então v é o gradiente de u no sentido da definição acima. 
Observe também que se u é diferenciável em algum aberto, então as componentes do gradiente 
segundo a definição acima coincidem com as respectivas derivadas parciais em quase todo ponto 
desse aberto. 

Denotamos por L 2 (Q) o conjunto {p : Q — > M m : pi G L 2 (Q),i = 1, . . . , m} com produto interno 
dado por 

ÍP, <?)l2 = / p- qdx (5.8) 



e norma ||p||l2 = (p,p) = J Q \p\ 2 dx = \\pi\\ 2 L 2 H h lb m ||^ 2 . 

Proposição 5.4 ~L 2 (Q) é um espaço de Hilbert. 

Prova: A verificação de que (J5.8)) define um produto interno é direta. Resta então mostrar que 
L 2 (Q) é completo. E claro que ||pí||l 2 < IIpIIl 2 para i = 1, . . . ,m e para todo p G L 2 . Seja (p n ) 
uma sequência de Cauchy em L 2 (Q). Temos que (pf) são sequências de Cauchy para i — 1, . . . , m 
e, como L 2 (Q) é completo, temos que p" —> pi G L 2 (Q). Tomando p = (pi, . . . ,p m ) G L 2 (Q), 

n— »oo 

temos que \\p n — 2 = \\Pi — Pi||r,2 H h ||p^ — Pm||r,2 — >• 0, logo p n — > p. Portanto, ~L 2 (Q) 

n— >oo n— >oo 

é completo. 

O espaço de Sobolev H l {Q) é definido por 

^{Q) = {u G L 2 {Q) : existe V« e este está em L 2 (Q)} . 
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Observe que faz sentido falar em gradiente de u na definição acima, pois L 2 (Q) C L l (Q) uma 
vez que Q é limitado. E claro que if 1 (<5) é um espaço vetorial e, sobre ele, podemos definir a 
forma bilinear 

(ui,u 2 ) = / (uiu 2 + V«i • Vw 2 ) dx. (5.9) 

E fácil verificar que essa forma bilinear define um produto escalar e que esse produto escalar gera 
a seguinte norma 

\\u\\ H i = y/(u,u) = ^\\u\\l+ \\u\\l (5.10) 
onde ||u||o = § Q u 2 dx = \\u\\ 2 L2 e = J Q \Vu\ 2 dx = ||Vu||^2- Mostraremos, a seguir, que 

Proposição 5.5 H X {Q), com produto interno dado por \5. fj\) . é um espaço de Hilbert. 

Prova: E suficiente mostrar que H 1 (Q), com a norma ||w||ííi = \/ (u, u), é completo pois, como 
já observamos, ()5.9|) define de fato um produto interno. Seja (w n ) n eN uma sequência de Cauchy 
em íT 1 (<5). Ponha v n = Vw n , n G N. Como ||w||z,2 < ||w||^i e ||Vu|| L 2 < jyi para toda 
u G i/ 1 (<5), temos que (u n ) e (v n ) são sequências de Cauchy em L 2 (Q) e L 2 (Q), respectivamente. 
Mas esses espaços são completos, logo u n — -> u G L 2 (Q) e v n —> v E L 2 (Q) quando n — > oo. 
Basta então mostrar que v = S/u, pois neste caso teremos u G H 1 ^) e, como \\u n — u\\ 2 HÍ = 
\\u n — u\\ 2 L2 + \\v n — v ||? 2 — ► 0, teremos também que it n w. Seja <p G Cq°(Q). Se = temos 

n— >oo n^oo 

(|5.7J) trivialmente; assuma então que 7^ 0. Seja e > 0. Faça M = ||0||_f/i > e tome n G N tal 
que || te — w n ||z,2 + \\v — v n \\- L 2 < e/M. Para cada i — 1, . . . , m, temos por Cauchy-Schwartz 



í (u-u n )^-(j)dx+ f(Vi-V? 
Jq o x í Jq 



)cf)dx 



< \\u - n n || L2 ||V0|| L 2 + \\v - v n \\tf\\(f>\\ L 2 < 

< (||« - M n || L 2 + || V - V n \\ h 2)M < 6. 



Agora note que, pela definição de gradiente, vale 

Q 9xí 



d 

u n -—ódx + I vfódx = 0. 



Logo, 

d 



u— — ódx + / Viódx 
dxi 



< e. 
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Como e > é arbitrário, temos que u e v satisfazem ()5.7|) e, portanto, v = Vtt. □ 

A desigualdade de Poincaré, que enunciamos abaixo 1 , será usada mais adiante no estudo da 
equação de Poisson com condição de contorno periódica. 

Teorema 5.6 (Desigualdade de Poincaré) Existe cq > 0, que depende apenas de Q C R m , 
tal que 

VueH\Q). 

Um campo de vetores v G L 2 (Q) é chamado de campo potencial, ou campo gradiente, se v = Vu 
para algum u G H 1 ^). 

Dado p G L 2 (Q), definimos o divergente de p como sendo o funcional linear dado por 

div p ■ (f) = (div p,cj)) = - / (p ■ V(j))dx, G ^(Q). 

Jq 

A desigualdade de Cauchy-Schwartz em L 2 (Q) implica que div p é contínua, pois 

(div p,(f>) < imi L 2||v0||l2 < iblMMIffi- 

Convém observar que, para funções diferenciáveis em um aberto de Q, as definições de divergente 
e campo potencial coincidem com aquelas definidas no cálculo vetorial. 

Outro exemplo de funcionais lineares e contínuos em H 1 (Q) são aqueles dados por u t— > fudx, 
onde / é um elemento de L 2 {Q). 

Espaços de funções periódicas 

Denotamos por □ o paralelepípedo em IR m , dado por □ = (0, l\) x • ■ ■ x (0,/ m ). Dizemos que 

uma função definida em IR m é periódica, de período □, se ela for uma função periódica em cada 

uma das variáveis x\, . . . , x m , e com períodos l\, . . . , l m , respectivamente. Neste caso, também 

dizemos que a função é □— periódica. 

1 A prova da Desigualdade de Poincaré pode ser encontrada em |12j . 



u dx < c 



udx 



iVwl dx 
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Se g(-) : M m — > R é uma função □— periódica e Lebesgue-integrável em □, definimos o valor 
médio de g por 

(g) = pj J^g(x)dx, 
onde |d| = I1I2 ■ ■ ■ l m é o volume do paralelepípedo □. 

A seguir, definiremos alguns espaços de funções que estarão relacionados, de alguma maneira, 
com as funções de M m em IR que sejam □— periódicas. Assim, o símbolo □ representa não apenas 
o domínio dado pelo hipercubo mas também que estamos tratando de funções periódicas. Em 
alguns casos, essa diferença vai além do domínio de definição das funções do espaço em questão, 
pois quando tratamos de espaços que levam em consideração a derivada das funções, estaremos 
impondo condições na fronteira deste domínio. 

Por exemplo, dado p G [l,+oo), denotamos por L P (D) às funções g(-) : M m — * E que sejam 
□—periódicas e tais que a norma 

\\9\\ P = (\g\ p ) 1/p 

seja finita. E claro que L 2 (\3) com produto interno (u, v ) = (uv) é um espaço de Hilbert. 
Neste caso, existe uma bijeção direta entre L 2 (Q) e L 2 {\3) com Q — □, ou seja, os espaços são 
essencialmente os mesmos, exceto porque as funções do primeiro espaço têm que ser estendidas 
de forma d-periódica para obtermos elementos do segundo. 

Denotamos por C°°(n) o conjunto das funções de M. m em M, D-periódicas e infinitamente diferen- 
ciáveis. Definimos como sendo o completamento de C°°(n) com respeito à norma || • ||#i. 
Em definimos o produto interno (u,v) = (uv + Vw • Vv) e norma 1 1 1 1 1 (□) = \/ (u, u). 

Observe que a diferença entre funções de H l (Q) e não reside simplesmente no domínio 

em que estão definidas, mesmo quando Q — □. De fato, se G C°°(]R m ) então claramente 
4>\q G C°°(Q). Em particular, se G C°°(D), vale 4>\q G C°°(Q). Entretanto, a recíproca 
é falsa: se estendemos uma função <p £ C°°(Q) para todo o W 71 de forma D-periódica, em 
geral não obtemos uma função de classe C°°, pois não temos qualquer informação sequer sobre a 
continuidade de 4> n a fronteira de □ - a função e suas derivadas teriam que "colar" nesta fronteira. 
O ponto aqui é que as funções de H 1 (Q) cujas extensões D-periódicas estão em H l (n) são somente 
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aquelas que podem ser arbitrariamente aproximadas em H 1 (Q) por funções de C°°(U). Podemos 
então pensar que é, em certo sentido, um subespaço fechado de 

De forma análoga ao que ocorre em H l (Q), temos que qualquer / £ L 2 {\3) define um funcional 
linear contínuo em ií 1 (D) dado por (/, u) = (fu) V u £ H l {\3). Da mesma forma, se g £ L 2 (D), 
então o funcional div g, dado por (div g,u) = — (g ■ Vu), é linear e contínuo em Uma 
consequência disso é que u i— » (u) é um funcional contínuo em eun (Vu) é uma aplicação 

contínua de em M. m . 

Definimos também o espaço 

Lj ot (D) = {/ : R m -> R m : 3 p £ M m , u £ ií" 1 (D) tais que / = p + Vu] . 

Claramente Lp 0t (d) é um espaço vetorial. Na verdade este espaço é também de Hilbert, mas não 
vamos utilizar este fato. 

Proposição 5.7 Todo campo de vetores v £ Lp 0t (D) é escrito como v = (v ) + Vu, u £ H l {\3). 

Prova: Primeiramente escrevemos v = p + Vu, p £ M. m ,u £ donde (v) = p + (Vu). O 

resultado desejado segue do fato de que (Vu) = V u £ if 1 (D). 



-ff 1 

Com efeito, se u £ if (□), temos que existe uma sequência (0 n ) em C°°(0) tal que <p n — > w, e, 

n— >oo 

como ii i— > (Vu) é contínua, temos que (V<f) n ) — > (Vu). 

Agora observe que, para qualquer função £ C°°(D), vale 
" k d 

——(j) dxi = (j)(xi, . . ., Xi-i, k, x i+ i . . . x m ) - <p(xi, . . ., Xi-i, 0, x i+í . . . x m ) = 

OXi 



e, pelo Teorema de Fubini aplicado à integral múltipla em temos que (^:> 

Repetindo esse argumento para i = 1, . . . , m, temos que (V(f>) = 0, completando a prova. 



□ 
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5.2.2 Lema de Lax-Milgram 



Seja V um espaço de Hilbert e a(u,v) : V x V — > IR um funcional bilinear. Suponha que a seja 
limitado, no sentido de que 

a(u, v ) < ^2||m||v|M|v (5-11) 
para todo u,v G V, onde v 2 > 0. Suponha também que a seja coercivo, isto é, 

a(it, «) > í/i||u||y, (5-12) 

para todo u E V , onde V\ > 0. 

Seja V a * o conjunto dos funcionais lineares e contínuos de V em R. Para cada u E V, a forma 
linear v t— > a(u,v) é contínua e, portanto, representa um elemento de V* , que denotaremos por 
Au. E claro que o operador A : V — >• V a * que leva m em Am é linear e satisfaz 

{Au, v) = a(u, v). 

onde (•, •) denota o produto escalar em V . Dada / G V* , considere a equação em u: 

a(u,v) = (f,v) \fveV. (5.13) 
isto é, dado / G V* , queremos resolver Au = f, com u G V . 

Teorema 5.8 (Lema de Lax-Milgram) O problema 115.1 fy) tem uma única solução u, que sa- 
tisfaz 

^ilMk<||/||v*<^IMIv. (5.14) 
Em outras palavras, uma forma bilinear limitada e coerciva define um isomorfismo entre V e V* . 

Prova: Seja / G AV C V* e seja u G V tal que Au = f. Por (|5.12|) temos 

i*i NI V < (Au,u) = (f,u) < ||/||y||«||y, 

donde segue a primeira desigualdade de (|5.14|) . Daí temos que A é injetivo, pois Au = Av implica 
A(u — v) = 0, logo ||« — v\\v < 0. Da injetividade de A segue a unicidade de soluções, quando as 
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mesmas existirem. Ainda supondo a existência de soluções, a segunda desigualdade em ()5.14|) é 
obtida de (|5.11|) pois, dado qualquer v G V, temos 

= \( Au i v )\ = \a(u,v)\ < v 2 \\u\\ v \\v\\ v . 



Resta mostrar que AV = V* . Das desigualdades acima segue que A : V — > AV C V* é um 
isomorfismo, e sequências de Cauchy em AV são imagens de sequências de Cauchy em V, que 
convergem porque V é completo. Portanto, AV é fechado em V*. 

Agora seja / G AV- 1 -. Pelo Lema de Riez, existe um único z G V tal que (/, v) = (z, v) V v G V. 
Temos que (Au, z) = (Au, f) = V u G V . Tomando u = z temos (Az, z) — e, por (|5.12[) . z — 0, 
logo / = 0. Como / G AV 1 - é arbitrário, temos que AV 1 - = {0}. Pelo Teorema da Projeção 
Ortogonal em Espaços de Hilbert, AV = V* mas, como AV é fechado, segue que AV = V*. 



5.2.3 O problema periódico 

Para tornar o texto mais limpo, vamos assumir, daqui em diante, que estaremos tratando da 
soma de 1 até m sempre que os índices i, j ou k aparecerem repetidos, a não ser que esteja 
explícito o contrário. 

Escrevemos A(x) = (a.ij(x) )™ =1 para denotar uma matriz cujos elementos : lR m — > IR são 
funções mensuráveis e limitadas. Fixado x G M m e dados Ç,rj G M m , definimos os produtos 
(A£)i = CLijÇj, (£A)j = Ci a ij- Escrevemos também rjAÇ = r)(A£) = (r)A)t; = rjiAij^j. 

Definimos 

\CA(x)r]\ 



\\A\\ = sup sup 

xeM m Ç,í76R m \{0} |ç|M 

de forma que 

\£A{x)ri\ < \\A\\ |eiM Vx.^er. 

Uma vez que os elementos de A satisfazem |ay(a;)| < M^, temos que ||^4|| < oo e Ylij & uma 
cota superior para 
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Daqui em diante, vamos assumir que A sempre satisfaz a condição de elipticidade 

para algum V\ > 0. Dizemos também neste caso que A é uniformemente elíptica. 
Considere a seguinte equação diferencial 

div(„4Vu) = -/o + div /, 

onde f eL 2 (D)e/e L 2 (D). 

Diremos que u G ií 1 (n) é uma solução fraca da equação acima se u satisfizer 

(V<P-AVu) = (f <f ) ) + (V<p-f), V^Gtf^D). (5.15) 

Mostraremos agora que o problema acima possui solução e a mesma é única no sentido de que 
duas soluções diferem por uma função constante, desde que (/o) = 0. 

Proposição 5.9 O problema \5.15\) possui solução u G se, e somente se, (fo) = 0. Neste 

caso, tal solução é única no sentido de que u — vé constante para qualquer outra solução v. Em 
particular, existe uma única solução satisfazendo (u) = 0. 

Prova: Suponha que existe uma solução u G H l {\3). Fazendo <p(x) = 1 V x G M m , temos V0 = 
e G if 1 (D). Substituindo em (|5.15j) . temos (fo) = 0. 

Reciprocamente, suponha que (/ ) = 0. Se u é uma solução então u + c é também solução para 
qualquer c G M, uma vez que somente Vw aparece em ()5.15|) . Assim, não é possível termos 
unicidade mais forte do que a do enunciado acima. Por outro lado, suponha que u, v são duas 
soluções do problema e escreva u = (u) + ú e v = (v) + v. Temos que ú e v são também soluções 
e satisfazem (ú) = (v) = 0. Como u — v — (u — v) + (ú — v), para que u — v seja constante basta 
mostrar que ú = v, ou seja, basta mostrar que existe uma única solução u satisfazendo (u) = 0. 

Para isso, defina V — {u G H 1 ^) : (u) = 0}, subespaço de iJ 1 (D). V é fechado porque u i— > (u) 
é contínua; logo V é um espaço de Hilbert. Considere o funcional bilinear a : V x V — ► IR dado 
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por a(u, 0) = (V0 • AVu). Vale 

a(u,u) = (Vu-AVu) > (z/i|Vu| 2 ) = v x -r=- í \Wu\ 2 dx, 

|Q| Jo 

onde z/i é a constante da condição de elipticidade de A. Da Desigualdade de Poincaré segue que 

cq / \Vu\ 2 dx > / u 2 dx, 



O JD 



uma vez que (u) = 0. Assim, temos 

a(u,u) > — (1 + c ) / |Vw| 2 cfe 

(l + c )|D| y n 



> 7 *\, , ( / \Vu\ 2 dx + / u 2 cfo 

" (l + co)|D| \Jn Jo 

Ví II 112 



(1 + co 

donde concluímos que a é coercivo. 

Sejam u,v E V. Temos 

\a(v,u) \ = | (Vu-AVv) | < ( \Vu-AVv\ ) < ||A|| ( |Vu| \Vv\ ) , 

donde 

|a(w,w)| 2 < ||A|| 2 (|Vw||Vt;|) 2 < ||A|| 2 (|V«| 2 )(|V*;| 2 ), 
logo |a(w,u)| < ||A|| 1 1 ix- 1| v\\v \\v- Portanto, a é limitado. 

Considere o funcional linear i— > (fo<f>) + (/ ■ V0), definido em V. Esse funcional é contínuo, pois 
sua primeira parte é dada por um elemento fo G L 2 {\3) e sua segunda parte é o divergente de 
— / G L 2 (D). Pelo Lema de Lax-Milgram, existe um único «gV tal que a(u, 0) = (V0 • AVu) = 
(fo4>) + (f ■ V0) para toda G V. Agora seja G ^(D). Temos = - (0) G V e V0 = V0. 
Logo, (V0 ■ AVu) = (V0 • AVu) = (j o 0) + (/ ■ V0) = ( f (<f> - (0)) ) + </ • V0) = (M) - 
{fo (0)) + (/ • V0) = (fo4>) + (/ • V0), a última igualdade valendo porque estamos assumindo 
que (/o) = 0. Portanto, (J5.15j) vale para toda G ií 1 (D). [— | 

5.2.4 Um problema periódico auxiliar 

Vamos considerar o problema 

div(Av) = 0, v G L 2 ot (D), (v) = A G W n . (5.16) 
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Todo elemento de Lp 0t (d) é escrito como v — X + Vit, u G ií 1 (n), sendo A = (v). Assim, 
obtemos o problema abaixo para u: 

à\v{A{X + Vu)) = 0, « 6 ^(D). 

Podemos reescrever o problema acima na forma (|5.15jl com / = —AX e fo = 0. Pela Proposição 15. 9| 
existe uma única solução u do problema acima, a menos de adição de uma constante. Como v 
depende apenas de Vw, adição de constante em u não altera a solução v, logo v é única. Em 
outras palavras, para cada A G M m , o problema (|5.16j) possui uma única solução v G Lp 0t (D). 

Por outro lado, v i— > div(„4.t>) é uma aplicação linear, logo seu núcleo é um subespaço vetorial. 
Assim, se v e w são duas soluções arbitrárias, isto é, div(Av) = div(Aw) = 0, e c G M é um 
escalar arbitrário, temos div(^4(^ + cu;)) = e (v + cw) = (u) + c (w) = X + c£, onde A = (v) e 
£ = (u>). Devido à unicidade da solução de (|5.16J) para cada parâmetro A fixo, podemos tomar o 
caminho inverso, ou seja, dados A, £ G M m e c G K arbitrários, temos que a solução do problema 
(|5.1fij) com parâmetro A + c£ é i; + cw. Assim, a solução it depende linearmente de A. 

Portanto, (Av) depende linearmente de A e, conseqúentemente, pode ser representada por 

(Au) = A°X, (5.17) 

onde A° é uma matriz real com elementos constantes. 

O problema periódico auxiliar e a determinação dos elementos da matriz A° em termos de A(x) 
é de importância central no estudo dos problemas homogeneizados - veja as seções subseqúentes. 

Vamos então tentar encontrar os elementos a^- de «4 o = (a° )™ =1 . Seja e 1 , . . . , e m a base canónica 
de M m . Para cada k — 1, . . . , m, a solução de ()5.16|) com A = e k é da forma v = e h + ViV^, G 
H\n). Temos 

a° = e* ■ A°e j = (é ■ Av) = {é ■ A{e ] + VN Ô )) = (a ió + e l ■ AVNj) , 

logo 

4 = K> + («^)- ( 5 - 18 ) 
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Por definição, iVjt é a solução de 



div(^(e fc + ViV fc )) = 0, N k e H\n) 

ou, no caso de A suave: 



(5 - l9) 



Proposição 5.10 A matriz Aq é elíptica. Além disso, se A(x) é simétrica para todo x, então 
A° é simétrica. 

Prova: Segue de (j5.16j) que 

V(f)±Av V^eif^D). 
Em particular, Av _L (v — A), logo (vAv) = (XAv). Disso segue que 

A ■ *4°À = A ■ (Av) = (A ■ Av) = (v ■ Av) > v x (|v| 2 ) > v x \ (v) \ 2 = u x \\\ 2 , 

ou seja, ^4° é elíptica. 

Para provar a segunda parte, consideramos o problema análogo 

divK4) = o, weL 2 pot (D), H=Cer, (5.20) 

cuja única solução w depende linearmente de £, e definimos a matriz constante C° por 

(wA) = ÇC°. (5.21) 
Do mesmo argumento de ortogonalidade, segue que 

Av±(w-Ç) e wA± (v - A), 

e, consequentemente, 

ÇA°\ = Ç (Av) = (ÇAv) = (wAv) = (wAX) = (wA) A = £C°A. 

A igualdade acima vale para todos A,£ G M m . Segue que A° = C°. Se definimos C° utilizando 
A T ao invés de A, obtemos (A°) T , implicando na seguinte propriedade: (*4°) T = (^4 T )°. Em 
particular, quando A(x) é simétrica para todo x, temos A° simétrica. 
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5.2.5 Homogeneização de equações parabólicas 



Nesta seção oferecemos uma prova construtiva para o Teorema 15. in utilizando o método conhecido 
como método de expansão assintótica. 

Este método é comum no estudo teórico de sistemas físicos e consiste em considerar formalmente 
a seguinte aproximação para u e : 

u e (x, t) th u°(x, t) + eui(x, y, t) + e 2 u 2 (x, y, t) + e 3 u 3 (x, y,t) H , 

onde y = e~ x x, as funções Ui, u 2 , u 3 , . . . são periódicas em y e não dependem de e e a função u 
depende apenas de x e t. 

A idéia do método é fazer com que o truncamento da expansão acima até o n-ésimo termo, 
que denotaremos por u e n , por um lado satisfaça "aproximadamente" a equação diferencial que 
estamos considerando e por outro lado seja bem aproximado pelo seu primeiro termo uq(x,Í), 
que não depende de e direta- nem indiretamente. 

Na prova que apresentamos abaixo vamos considerar a aproximação até segunda ordem em e, 
dada por 

u e 2 (x, t) = u (x, t) + eux(x, y, t) + e 2 u 2 (x, y, t), y = e~ x x. (5.22) 
Apresentamos abaixo o teorema mencionado e seu corolário. 

Teorema 5.11 Para cada e > 0, seja u £ a solução (clássica) do seguinte problema 

f p e u\ - divO^W) = 0, xe R m , t > 0, , , 

\ u e (x, 0) = cf>(x), 4> e C °°(R m ) [b ó) 

onde p e = p(e~ 1 x) e p(x) é uma função C 1 , positiva e periódica que satisfaz (p) = 1; A e = 
A(e~ 1 x) onde A(x) é uma matriz m x m periódica (de mesmo período de p), uniformemente 
elíptica, cujos elementos são funções C 1 de x. Então 

sup \u e {x,t) -u°(x,t)\ — ^0, Ví>0, (5.24) 

onde uq é a solução de 

u° t - div(^°VM°) = 0, u°(., 0) = (5.25) 
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e A° é a matriz constante dada por \5.1Tj) . 



Prova: A prova consiste em estudar a aproximação até segunda ordem para u € , dada por (|5.22|) . 

Primeiramente, definimos o operador 

d d d 

— = ^— + e — , z = l,...,m, 



e observamos que o operador div( v 4 e V) pode ser reescrito como 

Desenvolvendo o lado direito da igualdade acima, obtemos 

div(^ £ V) = e~ 2 A 1 + e~ 1 A 2 + A 3 , 

onde 



d 2 

<kj(y) 



A 1 


= ~( 

dyi \ 


A 2 


= A f 


A 3 


= a ió (y) 



9 2 



dxidxj 



Definindo o operador L e como 

^ = P e ^-div(^V), 
nós olhamos para a sua ação sobre w?>: 



- L e u\ 



/ du Q 

= e~\A lUl + A 2 u°) + í -p(y)_ + A lM2 + A 3 u° + A 2 u\ 

(ôu ôu \ 

A ãUl + A 2 u 2 - p(y) - ep(y) + eA 3 u 2 J . (5.26) 

Queremos determinar condições suficientes em u°, U\ e U2 para que L e u 2 seja arbitrariamente 
pequeno. 
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Começamos igualando o termo que multiplica e 1 a zero. Obtemos A\Ui + A2U = 0, ou seja 

d f d \ d ( . du°(x,t) 



Dado u°(x, t), vamos olhar para a igualdade acima como uma equação elíptica com incógnita u\ 
e variável independente y. Portanto, o problema acima é periódico em y, enquanto x e t são 
considerados como parâmetros. Pondo Nk(y), k = 1, . . . , m como as soluções de í[5.19|) . podemos 
tomar 

u 1 (x,y,t)=N k (y)^^- (5.27) 
ôx k 

como solução da equação acima. 

Agora vamos igualar a zero o termo que multiplica e°, obtendo assim um problema periódico 
para a incógnita u 2 e com variável independente y, dado por 

div y (V y u 2 {x,y,t)) = - l-p(y) — + A 3 {y)u {x, t) + A 2 (y) Ul (x, y, t) J . (5.28) 

Pela Proposição 15.91 a existência de uma solução u 2 para a equação acima depende da condição 

du° 
~dt 



du° 

-P(y) + A3U0 + A 2 u x ) = 0. (5.29) 



Mas 



A d 2 u°(x,t) 



A 3 u =aij(y)- 



dxidxj 



d ( du x \\ I d 2 u ± 
(A^) = (7- )) + ( aij- 



dyi V dxjj / \ dxídyj 
d 2 ui 

0+(a, 



v dxidyj 
d 3 N k u° 
13 dxidx k dy ó 
dNj\ d 2 u° 



dyk I dxidx 



.1 



onde na primeira igualdade substituímos a definição de A 2 ; na segunda, aplicamos o mesmo 
resultado contido na prova da Proposição 15.71 de que as derivadas de funções periódicas têm 
média zero; na terceira, substituímos ()5.27jl e, na quarta, apenas trocamos j por k. 
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Substituindo em (|5.29|) . temos 

du° _ I Õu° \ / dNA d 2 u° Õ 2 u° 



I du °\ ,a , \ I dN Á & u ° o 

{p-ãt) = {Á3Uo + Á2Ui) = \ a - + aik wj dífr s = a ^ 



dt \ dt I \ dyk I dxidxj %3 dxidxj 

onde a última igualdade vem de (|5.18|) . Portanto, u° deve satisfazer 

ua ' '0v7„.(h 



div(^ u Vu u ) = 



dt 

Ponha então u° como solução da equação acima com condição inicial u(-,0) = <p. Com u°, u\ e 
u 2 dados acima, acabamos de construir a aproximação (|5.22j) . 

Usamos a desigualdade triangular para obter 

sup \u e — u \ < sup \u e — u e 2 \ + sup \u e 2 — u°\, 

R m x[0,T] R m x[0,T] R m x[0,T] 

e vamos tentar controlar cada termo separadamente. 

O primeiro termo depende essencialmente de \L t u\(x,t)\ e de \u e (x, 0) — w|(a;, 0)|. Isso vem do 
Corolário do Apêndice lfí.2l aplicado à função u e — u\ e do fato de que u £ , sendo solução, 
satisfaz L e u e = 0: 

sup \u e — u\\ < sup \u e (x, 0) — u 2 (x, 0)| + (minp)~ 1 T sup \L t u 2 (x,t)\. (5.30) 

M m x[0,T] M m R m x[0,T] 

O teorema está provado se mostramos que, dado T > 0, valem os limites abaixo, quando e — >• 0: 

sup |u|-u°|->0 (5.31) 

R m x [0,T] 

sup|w e (x,0) -u e 2 (x,0)\ ^0 (5.32) 
sup \L e u 2 \ = sup |£ e (tÍ2 — u e )\ —> 0. (5.33) 

M m x[0,T] R m x[0,T] 

Comecemos por (|5.31|) . E claro que u° e suas derivadas são limitados em M." 1 x [0, T]. Consideran- 
do que as funções Nk(y) são periódicas e de classe C 1 , portanto limitadas, temos que u±, dado por 
(J5.27|) . e u 2 , dado pela única solução de (|5.28|) satisfazendo (u 2 ) = 0, são limitados em IR m x [0, T}. 
Segue então que u 2 , dada por ()5.22j) . converge uniformemente em W m x [0,T] para u° quando 
e — > 0. O limite (|5.32|) é uma conseqúência direta de (|5.31|) pois u e (x, 0) = u°(x, 0) = 0(x). Para 
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obtermos ()5.33|) . primeiro note que L e u\ é dado apenas pelo termo ()5.26|) . uma vez que u°, U\ e 
u 2 foram definidos de forma a anular os demais termos. Mas da mesma forma que U\ e u 2 são 
limitados em M. m x [0, T], também o são suas derivadas com respeito a x, y e t, de forma que 
L £ u 2 converge uniformemente para zero em IR m x [0,T]. Observe que estamos nos baseando no 
fato de que a solução u 2 da equação elíptica ()5.28jh assim como suas primeiras derivadas, são 
uniformemente limitadas, uma vez que Uq, sendo solução da equação do calor, também tem suas 
derivadas limitadas e u\ é dado por 1)5.27(1 . 

Os limites acima, juntamente com a estimativa ()5.30() . resultam em 

lim sup \u e — u°\ =0. 

e ^°R m x[0,T] 

Agora note que ()5.24j) é consequência direta deste último limite e a prova está concluída. 
O corolário abaixo justifica o argumento formal utilizado no Capítulo HJ 

Corolário 5.12 Para cada e > 0, seja u e a solução clássica do seguinte problema 

f u\ = [1 + ^(e^aOfe, x e R, t > 0, 
\ u e (x, 0) = <j>(x), (j) G C^(R m ) 

onde g é suave e periódica e \x satisfaz [1 + fig(y)] > V y G IR. Seja u° solução de 

u° t =au° xx , M °(-,í = O) = 0, (5.34) 
onde cr = ([1 + fig^y)]^ 1 ) 1 . Então 

sup \u € {x,t) -u°(x,t)\ — ^0, Ví>0. (5.35) 



Prova: Podemos reescrever a equação para u e na forma ()5.23|) com 

p(y) = \^ + w(y)Y 1 

e A(x) constante igual a 1. Pelo Teorema 15.111 vale o limite ()5.35|) sendo u° a solução de 

(p)u° t =u° xx , M °(-,í = O) = 0. 

Mas, pela definição de a, temos (p) = a^ 1 , portanto u° é também a solução de 1)5.34)1 e não há 
mais nada a provar. |— I 
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Capítulo 6 
Apêndices 



6.1 A transformada de Fourier 

Neste apêndice vamos definir a transformada de Fourier e enunciar as suas propriedades que são 
utilizadas no Capítulo |2] As demonstrações dos resultados abaixo podem ser encontrados no 
livro de Iório [TH] • Representamos por dx a medida de Lebesgue na reta. 

Definição 6.1 (Transformada de Fourier) Seja f um elemento de L 1 (R). Definimos a trans- 
formada de Fourier de f através da fórmula 

Hf{x)) = í f(x)e~ íkx dx = />). (6.1) 
Jm 

atua linearmente em L X (]R) e além disso, para qualquer elemento / G L 1 (M) e a G M, 

F{f{ax)} = ±-J{k/a). (6.2) 
\a\ 

Outra propriedade é que se / G L 1 (1R), então / é contínua e limitada. 
Para funções dentro da classe de Schwartz iS(M) 1 , vale o seguinte resultado: 

Teorema 6.2 Se f E S então f G S e 

f'(k) = -ikf(k). (6.3) 

1 S(R) = {/ : R -> C/\x m D n f(x)\ < M(m,n) V m, n e N}, onde D n f(x) denota a n-ésima derivada de / no 
ponto x. 
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Além disso, para quaisquer f,g £ S: 

{f,g)=2*{f,g), (6.4) 
onde (•, •) denota o produto interno de L 2 (R). Em particular, para todo f £ S, 

^ = ^11/1112. (6.5) 



Observação: Como S é um subespaço denso em L 2 (R), a definição de J 7 , pode ser estendida ao 
espaço de Hilbert L 2 (R) ainda satisfazendo a relação (|fi.5j) . A transformada inversa de Fourier 
está bem definida em L 2 (R) e a mesma é dada pela fórmula de inversão abaixo. 

Teorema 6.3 (Transformada inversa de Fourier) Para todo f £ L 2 (R) ílL^l), temos que 
f G L 2 (M) e, além disso: 

1 r R « 

f(x) = — lim / f(k)e ikx dk, (6.6) 



2tt r- 



R 



o limite valendo em L 2 (R). Para uma subseqúência R n — >• oo o limite também vale q.t.p. 

Definição 6.4 (Convolução) Sejam f,g £ L 1 (R). Definimos a convolução de f com g pela 
integral 

(f*9)(x)= / f(x -y)g(y)dy. (6.7) 



A convolução de duas funções / e g satisfaz às seguintes propriedades: f*g = g*fe (/ * g) * h = 
f * (g * h). O teorema da convolução nos diz que: 

Teorema 6.5 (Teorema de convolução) Para todo f,g £ <S(R): 

F{f*g}(k) = f{k)g{k), (6-8) 

ou seja, 

F~ 1 {fg}(x) = (f*g)(x). (6.9) 

Para finalizar, lembremos que a distribuição gaussiana se transforma da seguinte maneira: 



^je^j = VÃo^e- a k , (6.10) 
onde a é uma constante estritamente positiva. 
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6.2 O princípio do máximo 



Neste apêndice vamos estudar propriedades sobre o máximo de funções que satisfazem uma 
equação (ou melhor, uma inequação) diferencial dada por um operador parabólico, que consiste 
em tomarmos o divergente do produto de uma matriz uniformemente elíptica 2 pelo gradiente 
da função, subtraído de sua derivada no tempo. Trata-se de uma das várias formulações do 
conhecido Princípio do Máximo para equações diferenciais parabólicas. 

Vamos enunciar o teorema abaixo sem no entanto demonstrá-lo 3 . 

Teorema 6.6 Seja u : M, m x [0, oo) — > E de classe C 2 satisfazendo a seguinte desigualdade 

(ô 2 u \ ( ôu \ õu 

em Et = M. m x (0, T]. Suponha que os coeficientes de L são contínuos e limitados e que a matriz 
(aij) é uniformemente elíptica. 

Nessas condições, se a função u é limitada em Et e u(x, 0) > V x E M. m então u(x,t) > 
V (x,t) e E T . 

Corolário 6.7 Seja p : R m — > M limitada de classe C 1 tal que inf K m p(x) > 0. Então o Teore- 
ma \6.b\ continua válido se definimos 

LN = |-(a í ,(x,t)^-p(x)|, (6.11) 

desde que os elementos da matriz (aij) sejam de classe C 1 e suas derivadas sejam limitadas. 

Prova: Desenvolvendo a primeira parte de L obtemos 

Definimos L[u](x,t) = p(x) _1 L[m](x, t), cujos coeficientes ainda satisfazem as hipóteses do Teo- 
rema EHU O corolário segue então do fato de que L[u] < ^ L[u] < 0. |— | 
2 Veja a condição de elipticidade definida na págma|HHl 

3 Uma versão mais geral desse teorema pode ser encontrada em p. 43, Thcorem 9. 
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Corolário 6.8 Seja u : M. m x [0, oo) — > M. uma função de classe C 2 tal que L[u] > em Et = 
K' m x (0, T] e u é limitada em Et, onde L[-] é definido por H6.ll]) . Então 

sup u(x,t) < supu(x, 0). 

R m x[0,T] R m 

Prova: Defina v(x, t) = [sup E m u(x, Q)]—u(x, t). Claramente v é de classe C 2 e satisfaz v(x, 0) > 
para todo x G M m . Além disso, L[v] = —L[u] < e v é limitada em Et- Logo, v(x, t) > em 
Et- Substituindo, obtemos que u(x,t) < [sup Rm u(x, 0)] para todo par (x,t) G Et- □ 



Corolário 6.9 Seja v : R m x [0, oo) 
JOTI) . Então 



duas vezes diferenciável e considere L[-] definido por 



Prova: Defina 



Segue que 



sup \v(x, í)| < sup \v(x, 0)| + (inf p) 1 T sup \Lv(x, t)\. 

m x[0,T] K m Rm M m x[0,T] 



u(x, t) = — ( (inf p) 1 sup \Lv{x, t)\ ] í + v. 

R m x [0,T] 



(infp) 1 sup \Lv(x,t)\ ]t 

R m x[0,T] 



Lv 



P 



> 



sup \Lv(x,t) \ + Lv 

inf p R m x [ T ] 

p 

, sup \Lv(x, t)\ — \Lv\ 

mf P Rm x [Q,T] 



> sup |Lu(x, í)| — |Lu| 

R m x [0,T] 

> o. 

Do corolário anterior segue que 

v (x, t) < sup \v(x, 0)| + (infp) -1 sup |Lf(a;,í)| }í 
M m y R m x[o,T] y 

para todo t G [0, T]. Aplicamos para — v o resultado que acabamos de obter e concluímos que 



\v(x, t)\ < sup \v(x, 0)| + I (inf p) 1 sup |Lt>(x, t)\ \ t 

R m y R m x[0,T] 

para todo t G [0, T]. A estimativa desejada segue direto desta última. 



□ 



108 



Referências Bibliográficas 



[1] W. F. Ames. Numerical Methods for Partial Differential Equations. Academic Press, 
3 a edição, 1992. 

[2] A. Bensousan, J. Lions, G. Papanicolaou. Asymptotic Analysis of Periodic Structure. North 
Holland, Amsterdam, 1978. 

[3] G. A. Braga, F. Furtado, J. M. Moreira, L. T. Rolla. Renormalization group analysis of non- 
linear diffusion equations with time dependent diffusion coemcients. Artigo em preparação. 

[4] G. A. Braga, F. Furtado, J. M. Moreira, L. T. Rolla. Renormalization group analysis of 
nonlinear diffusion equations with periodic coemcients. Multiscale Modeling and Simulation, 
1 (4), pp. 630-644, 2003. 

[5] G. A. Braga, J. M. Moreira, L. T. Rolla. Análise assintótica de soluções de EDP's via grupos 
de renormalização. Em Atas do 57° Seminário Brasileiro de Análise, pp. 213-280, 2003. 

[6] J. Bricmont, A. Kupiainen. Renormalizing partial differential equations. Constructive 
Physics, 446, pp. 83-115, 1995. 

[7] J. Bricmont, A. Kupiainen, G. Lin. Renormalization group and asymptotics of solutions 
of nonlinear parabolic equations. Communications in Pure and Applied Mathematics, 47, 
pp. 893-922, 1994. 

[8] L. Chen, N. Goldenfeld. Numerical renormalization group calculations for similarity solutions 
and travelling waves. Physical Review E, 51, pp. 5577-5581, 1995. 

[9] G. Dagan. Flow and Transport in Porous Formations. Springer- Verlag, New York, 1989. 



109 



[10] C. R. de Oliveira. Introdução à análise funcional. Publicações Matemáticas. IMPA, Rio de 
Janeiro, 2001. 

[11] G. Duro, E. Zuazua. Large time behavior for convection-diffusion equations in R n with 
periodic coefficients. Journal of Differential Equations, 167, pp. 275-315, 2000. 

[12] L. C. Evans. Partial Differential Equations, volume 19 de Graduate Studies in Mathematics. 
AMS, 1998. 

[13] A. Friedman. Partial Differential Equations of Parabolic Type. Krieger, 1983. 

[14] F. Furtado, J. Glimm, B. Lindquist, F. Pereira. Multi-length scale calculations of mixing 
length growth in tracer floods. Em F. Kovarik, editor, Proc. of the Emerging Technolo- 
gies Conference, pp. 251-259. Institute for Improved Oil Recovery, University of Houston, 
Houston TX, 1990. 

[15] J. Glimm, B. Lindquist, F. Pereira, Q. Zhang. A theory of macrodispersion for the scale up 
problem. Transport in Porous Media, 13, pp. 97-122, 1993. 

[16] N. Goldenfeld. Lectures on Phase Transitions and the Renormalization Group. Addison- 
Wesley, Reading, 1992. 

[17] D. Halliday, R. Resnick, K. S. Krane. Physics. John Wiley, 5 a edição, 2000. 

[18] R. Iório Jr., V. M. Iório. Equações Diferenciais Parciais: Uma Introdução. IMPA, Rio de 
Janeiro, 1988. 

[19] V. Isaia. Intermediate Asymtotic Behavior of Nonlinear Parabolic PDEs via a Renormaliza- 
tion Group Approach: a Numerical Study. Tese de doutorado, Department of Mathematics, 
University of Wyoming, Laramie, Wyoming, 2002. 

[20] V. V. Jikov, S. M. Kozlov, O. A. Oleinik. Homogenization of Differential Operators and 
Integral Functionals. Springer Verlag, Berlin, 1994. 

[21] R. J. LeVeque. Numerical Methods for Conservation Laws. Lectures in Mathematics ETH 
Zurich. Birkháuser Verlag, Boston, 2 a edição, 1992. 



110 



[22] J. M. Moreira. O Comportamento Assintótico de Soluções da Equação do Calor não-Linear 
Via Grupos de Renormalização. Tese de mestrado, UFMG, Belo Horizonte, Minas Gerais, 
2002. 

[23] A. Naddaf, T. Spencer. On homogenization and scaling limit of some gradient perturbations 
of a massless free field. Comm. Math. Phys., 183, pp. 55-84, 1997. 

[24] G. Papanicolaou. Diffusion in random media. Surveys in Applied Mathematics, pp. 205-255, 
1995. 

[25] D. W. Peaceman. Fundamentals of Numerical Reservoir Simulation. Elsevier, New York, 
1977. 

[26] L. Tartar. Homogénéisation. Cours Pecout ou College de France, 1977. 



111 



Errata 



Esta errata contém os erros encontrados até 08/05/2006. 

Pede-se comunicar quaisquer outros erros matemáticos, tipográficos ou gramaticais através do 
endereço eletrônico leorolla@impa.br. 



Página IH01 Na equação (j4.11j) deve-se ler 



ln Ar, ~ ln 



A ' 
(lnL)«_ 



— ahxn. 
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